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Abstract: Using holography, we study the evolution of a spatially homogeneous, far 
from equilibrium, strongly coupled A/" = 4 supersymmetric Yang-Mills plasma with a 
non-zero charge density or a background magnetic held. This gauge theory problem 
corresponds, in the dual gravity description, to an initial value problem in Einstein- 
Maxwell theory with homogeneous but anisotropic initial conditions. We explore the 
dependence of the equilibration process on different aspects of the initial departure 
from equilibrium and, while controlling for these dependencies, examine how the equi¬ 
libration dynamics are affected by the presence of a non-vanishing charge density or an 
external magnetic held. The equilibration dynamics are remarkably insensitive to the 
addition of even large chemical potentials or magnetic helds; the equilibration time is 
set primarily by the form of the initial departure from equilibrium. For initial devia¬ 
tions from equilibrium which are well localized in scale, we formulate a simple model 
for equilibration times which agrees quite well with our results. 
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1 Introduction 

The discovery of gauge/gravity duality (or “holography”) has enabled the study of 
previously intractable problems involving the dynamics of strongly coupled gauge the- 
oriesd In the limit of large gauge group rank Nc, and large ‘t Hooft coupling A, the 
strongly coupled quantum dynamics of certain gauge held theories may be mapped, 
precisely, into classical gravitational dynamics of higher dimensional asymptotically 

^See, for examples, refs. [1-3] and references therein. 
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anti-de Sitter (AdS) spacetimes [4-6]. Numerical studies of the resulting gravitational 
dynamics can shed light on poorly understood aspects of the quantum dynamics of 
strongly coupled gauge theories. 

Using the simplest example of gauge/gravity duality, applicable to maximally su¬ 
persymmetric SU{Nc) Yang-Mills theory [Af = 4 SYM), this approach has been applied 
to a succession of problems of increasing complexity involving far from equilibrium dy¬ 
namics. These include homogeneous isotropization [7-9], colliding shock waves [10-15], 
and turbulence in two-dimensional fluids [16, 17]. A detailed presentation of the meth¬ 
ods used in most of these works is available [18]. 

In this paper, we extend previous work on the dynamics of homogeneous but 
anisotropic Af = 4 SYM plasma [7-9]. We examine the influence on the equilibration 
dynamics of a non-zero global U{1) charge density, or a background magnetic held. 
Inclusion of these effects is motivated by the physics of relativistic heavy ion colli¬ 
sions [19-21]. Hydrodynamic modeling of near-central events clearly indicates that the 
baryon chemical potential fis in the mid-rapidity region is signihcantly smaller than 
the temperature, but not by an enormously large factor at RHIC energies.^ Hence, it 
is desirable to understand the sensitivity of the plasma equilibration dynamics to the 
presence of a baryon chemical potential and associated non-zero baryon charge density. 
Similarly, it is clear that large, but transient, electromagnetic helds are generated in 
heavy ion collisions. A growing body of work [24-28] suggests that electromagnetic 
effects may play a signihcant role despite the small value of the hne structure con¬ 
stant. Electromagnetic effects on equilibrium QCD properties are also under study 
using lattice gauge theory [29-32]. 

The large Nc, strongly coupled Af = 4 SYM plasma we study is, of course, only a 
caricature of a real quark-gluon plasma. But it is a highly instructive caricature which 
correctly reproduces many qualitative features of QCD plasma (such as Debye screen¬ 
ing, hnite static correlation lengths, and long distance, low frequency dynamics de¬ 
scribed by neutral fluid hydrodynamics). Moreover, in the temperature range relevant 
for heavy ion collisions, quantitative comparisons of bulk thermodynamics, screening 
lengths, shear viscosity, and other observables show greater similarity between Af = 4 
SYM and QCD than one might reasonably have expected [33, 34]. Since the composi¬ 
tion of a plasma depends on the chemical potentials, or associated charge densities, of 
its constituents, studying the dependence of the equilibration dynamics on a conserved 
charge density provides a simple means to probe the sensitivity of the dynamics to the 
precise composition of the non-Abelian plasma. This, in small measure, may help one 

^Inferred values of ^j-b/T at chemical freeze-out are about 0.15 for RHIC collisions at -^snn = 200 
GeV, and roughly 0.005 for LHC heavy ion collisions with yY/Uv = 2.8 TeV [22, 23]. 
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gauge the degree to which A/" = 4 SYM plasma properties can be extrapolated to real 
QCD plasma. At the very least, strongly coupled A/' = 4 SYM theory provides a highly 
instructive toy model in which one may explore, quantitatively, non-trivial aspects of 
non-equilibrium gauge held dynamics.^ 

The remainder of the paper is organized as follows. Section 2 summarizes necessary 
background material. This includes the coupling of an Abelian background gauge 
held to a U{1) subgroup of the SU{4:)ji global symmetry group of A/" = 4 SYM. This 
17(1) symmetry may be regarded as analogous to either the baryon number U{1)b or 
electromagnetic 17(1)em havor symmetries of QCD. Turning on a background magnetic 
held implies an enlargement of the theory under consideration from A/" = 4 SYM to 
AT = 4 SYM coupled to electromagnetism (which we abbreviate as SYM-I-EM). The 
combined theory is no longer scale invariant; this has important implications which 
we discuss. This section describes the 5D Einstein-Maxwell theory which provides the 
holographic description of the states of interest, presents our coordinate ansatz (based 
on a null slicing of the geometry), and summarizes relevant portions of the holographic 
dictionary relating gravitational and dual held theory quantities. This section also 
records the reduced held equations which emerge from our symmetry specializations, 
describes the relevant near-boundary asymptotic behavior, and summarizes properties 
of the static equilibrium geometries to which our time dependent solutions asymptote 
at late times. 

The following section 3 briehy describes our numerical methods, which are based 
on the strategy presented in ref. [18]. When studying states with a non-zero charge 
density (but no background magnetic held) appropriate numerical methods for asymp¬ 
totically AdS Einstein-Maxwell theory are immediate generalizations of methods which 
have previously been found to work well for pure gravity. However, the inclusion of 
a background magnetic held induces a trace anomaly in the dual quantum held the¬ 
ory which, in the gravitational description, manifests in the appearance of logarithmic 
terms in the near-boundary behavior of helds. Such non-analytic terms degrade the 
performance of spectral methods, on which we rely, and necessitate careful attention 
to numerical issues. Section 3 also describes the specihcs of our chosen initial data. 

Results are presented in section 4. We focus on the evolution of the expectation 
value of the stress-energy tensor. We hrst discuss the sensitivity of the equilibration 
dynamics to features in the initial data and, in particular, examine the extent to which 
the evolution shows nonlinear dependence on the initial departure from equilibrium. 
We hnd that only disturbances in the geometry originating deep in the bulk, very 

^Previous work examining thermalization in plasmas with non-zero chemical potential (not involv¬ 
ing numerical solutions of far from equilibrium geometries) includes refs. [35-38]. 
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close to the horizon, generate signihcant nonlinearities. This is broadly consistent with 
earlier work [8, 9]. However, for a very wide variety of initial disturbances, including 
ones which generate extremely large pressure anisotropies, we hnd remarkably little 
nonlinearity in the equilibration dynamics, often below the part-per-mille level. 

We then present comparisons of the equilibration dynamics as a function of the 
charge density or background magnetic held. We focus on comparisons in which the 
form of the initial departure form equilibrium and the energy density, or the equilib¬ 
rium temperature, is held hxed while either the charge density or magnetic held is 
varied. These comparisons reveal surprisingly little sensitivity to the charge density, or 
magnetic held, even at early times when the departure from equilibrium is large. 

We verify the late time approach to the expected equilibrium states, and extract 
the leading quasinormal mode (QNM) frequency from the late time relaxation. Quasi¬ 
normal mode frequencies extracted from our full nonlinear dynamics are compared, 
where possible, with independent calculations of QNM frequencies based on a linearized 
analysis around the equilibrium geometry. This provides a useful check on our numerical 
accuracy. 

We dehne an approximate equilibration time based on the relative deviation of the 
pressure anisotropy from its equilibrium value, and examine the dependence of this 
time on charge density or external magnetic held. Once again, changes in this quantity 
are largest for initial disturbances which originate very close to the horizon, but the 
overall sensitivity of the equilibration time to the charge density or magnetic held is 
remarkably modest. 

The hnal section 5 discusses and attempts to synthesize the implications of our 
results. We present a simple model of equilibration times, for initial disturbances 
which are well localized in scale, which agrees rather well with our numerical results 
(but becomes less accurate for disturbances localized extremely close to the horizon). 
We end with a few concluding remarks.'^ 

2 Ingredients 

2.1 A/^ = 4 SYM in an external field 

We study maximally supersymmetric SU{Nc) Yang-Mills theory (A/" = 4 SYM) on four 
dimensional Minkowski space when the conserved current for a U{1) subgroup of the 

^As this paper neared completion, we learned of the somewhat related work by A. Buchel, M. 
Heller, and R. Myers [39]. These authors examine quasinormal mode frequencies in W = 2* SYM 
and argue that, in this non-conformal deformation of W = 4 SYM, the longest equilibration times are 
largely set by the temperature with little sensitivity to other scales. 
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S'f/(4)/j global symmetry group either (a) has a non-vanishing charge density, or (b) 
is coupled to a background Abelian gauge held describing a uniform magnetic held. 
The embedding of the U{1) symmetry is chosen such that the U{1) commutes with an 
SU{3) subgroup of the SU{4:)r global symmetry. 

The coupling to the external held has the usual form^ 

^ = 5 sym + j d^xr{x)A^^\x ), ( 2 . 1 ) 

where j°‘{x) is the conserved U{1) current normalized such that the four Weyl fermions 
of A/" = 4 SYM have charges {-|-3, -1,-1, — 1}/a/3 and the three complex scalars have 
charge -1-2/a/3. The overall factor of 1/a/3 in these charge assignments has no physical 
signihcance, but is chosen so that the trace anomaly and electromagnetic beta function 
(induced when this current is gauged) have convenient coefficients, as will be seen 
below.® The background U{1) gauge held A®^*(x) we take to have the form 

+ ( 2 . 2 ) 

with /i the chemical potential which, in equilibrium, will be conjugate to the charge 
density j®, and B the amplitude of a constant magnetic held pointing in the x^ direction. 
Although it should be straightforward to study dynamics when both the charge density 
j® and magnetic held B are non-zero, in this paper we focus for simplicity on the cases 
of either a non-zero charge density with vanishing magnetic held, j® / 0 and B = 0, or 
non-zero magnetic held with vanishing charge density, = 0 and B ^ 0. 

With a non-zero magnetic held B in the x^ direction, changes in the background 
gauge held under a translation in the x^ or x"^ directions, or a rotation in the x^- 
plane, can be compensated by a suitable U{1) gauge transformation. Hence, the 
theory retains full spatial translation invariance as well as rotation invariance in the 
x^-x‘^ plane. 

We will be interested in initial states which: (i) have non-trivial expectation values 
(r"^(x)) and (j"(x)) for the stress-energy tensor and U{1) current density, respectively; 
{ii) are invariant under spatial translations as well as 0(2) rotations in the x^—x'^ plane; 
and {in) are invariant under the SU{3)r subgroup of the SU (4)/^ global symmetry which 
commutes with our chosen f/(l). 

Since all A/” = 4 SYM helds transform in the adjoint representation of the SU{Nc) 
gauge group, the stress-energy and 0(1) current expectation values both scale as 0{N^) 

®We use a mostly-plus Minkowski space metric, = diag(—1,-1-1,-tl,-tl). 

®These charge assignments are 1 /a/ 3 times those used in ref. [40]. Overall rescaling of these charge 
assignments has implications for the holographic description which are noted below in footnote 11. 
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in the large limit. For later convenience, we define a rescaled energy density e and 
charge density p, via 



O 

III 

(2.3) 

K ^ {N^- 

-l)/(2vr^) . 

(2.4) 


A/” = 4 SYM is a conformal field theory with a traceless stress-energy tensor. Adding 
a chemical potential p introduces a physical scale, but does not modify the microscopic 
dynamics of the theory and hence does not affect the tracelessness of the stress-energy 
tensor. In contrast, introducing an external magnetic field does affect the microscopic 
dynamics and, in particular, generates a non-zero trace anomaly,' 

= = (2.6) 

The trace anomaly generated by the external magnetic field implies that the the¬ 
ory is no longer scale invariant. For example, the ground state energy density, as a 
function of magnetic field, need not have the simple form of some pure number times 
This will be seen explicitly below. The trace anomaly implies that there must be 
logarithmic dependence on a renormalization point. To interpret this dependence, it 
is appropriate to adopt the perspective that adding an external magnetic field means 
that the theory under consideration has been enlarged — it is now TV = 4 SYM coupled 
to U{1) electromagnetism (SYM-I-EM). The complete action of the theory is the SYM 
action, minimally coupled to the U{1) gauge field, plus the Maxwell action for U{1) 
gauge field. 


‘S’sYM+EM — coupled + S'em , (2.6) 

with 

Sem^-I d^x^Fl,. (2.7) 

The electromagnetic coupling (having been scaled out of covariant derivatives) ap¬ 
pears as an inverse prefactor of the Maxwell action. We regard the electromagnetic 
coupling as arbitrarily weak. Hence, quantum fluctuations in the U{1) gauge field 

^We define the external gauge field such that no factor of an electromagnetic gauge coupling appears 
in the interaction (2.1), in our 17(1) covariant derivatives, or in the trace anomaly (2.5). The coefficient 
— in the trace anomaly (2.5) equals the EM beta function coefficient b^, given below in eq. (2.9). 
(Note that the sign of the trace anomaly depends on the metric convention in use.) 
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are negligibly small, allowing us to view the EM gauge field as a classical background 
field.® 

However, just as in QED, fiuctuations in the SYM fields which are electromag- 
netically charged will cause the electromagnetic coupling to run with scale. The 
associated renormalization group (RG) equation for the inverse coupling has the usual 
form, 

e“^ = /?i/e 2 (e"^) = -bo + O(e^), (2.8) 

with the one-loop beta function coefficient® 

bo = K^l + T 

a a 

Here, g" = (3, —1, —1, —1 )/a/ 3 and = (2, 2, 2)/\/3 are the charge assignments of the 
four Weyl fermions and three complex scalars, respectively. Integrating this renormal¬ 
ization group equation leads, as usual, to 

l/e^(/i) = 6oln(AEM/fi) 0[ln(ln Aem//^)] , (2.10) 

with the RG invariant scale Aem denoting the Landau pole scale where the (one loop 
approximation to the) electromagnetic coupling diverges. 

The total stress-energy tensor derived from the combined action (2.6) will equal the 
AT = 4 SYM stress-energy tensor, augmented with minimal coupling terms to the EM 
gauge field, plus the classical Maxwell stress-energy. An essential point, however, is that 
while the total stress-energy tensor is well-defined, partitioning the stress-energy tensor 
into separate SYM and EM contributions is inherently ambiguous, as the individual 
pieces depend on the renormalization point. We define 

, ( 2 . 11 ) 

®In an arbitrary background SU{4) gauge field, the divergence of the 5'{7(4 )h current acquires an 
anomalous contribution, oc . This anomaly, when specialized to our chosen U{1) 

subgroup, is proportional to the sum of the cubes of our fermion charges and is non-zero, X)a(9f )^ = 
8 /a/3. To make the combined SYM-|-EM theory well defined, one could add to the theory additional 
fermions, charged under the U{1) but with no SYM interactions, which would cancel this U{1) anomaly. 
As we are not concerned with quantum fiuctuations in the C/(l) gauge field, the presence of this U{1) 
anomaly (in the absence of compensating spectators) is irrelevant for our purposes. 

non-renormalization theorem in supersymmetric Af = 4 SYM implies that the short distance 
behavior of the current-current correlation cannot depend on the ‘t Hooft coupling A [41]. This 
implies that the leading EM beta function coefficient bg does not depend on A, and hence may easily 
be evaluated in the A —>■ 0 limit. 


= K. 


(2,9) 
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with^o 


( 2 . 12 ) 


and 

= ^SYM,min. coupled(A'') • (2.13) 

The partitioning (2.11) of the stress-energy tensor puts all quantum corrections other 
than the running of the EM coupling into the SYM contribution AT^Yj^(/i). The scale 
dependence must, of course, cancel between the two terms because the total stress- 
energy tensor is a physical quantity. Therefore, the scale dependence in the SYM 
contribution to the stress-energy must simply compensate the known running of the 
inverse electromagnetic coupling (2.8) in the Maxwell stress-energy tensor (2.12), 

T“S(a) = bo [F“-'F^ - ybF>‘-F,,] . (2.14) 

Specializing to zero temperature states in a constant static magnetic held jB, the 
scale dependence (2.14) plus dimensional analysis implies that the SYM contribution 
to the ground state energy density is a non-analytic function of magnetic held, 

£(fi) = CoB" - il‘) = lBMn[B-(f,)/|B|], (2.15) 

with Co some pure number. (Here and henceforth, when considering physics in a non¬ 
zero magnetic held c(/r) = AT^ymM /^ denotes the SYM portion of the rescaled energy 
density.) In the second form of eq. (2.15), the analytic term has been absorbed by 
dehning a scale dependent “hducial” magnetic held amplitude, 

H*(/i)=/iV"Y (2.16) 

Note that the ground state energy acquires a simple quadratic form when the renor¬ 
malization point is chosen to scale with the magnetic held, = CqH^. Our 

numerically determined value for the coefficient cq is given below in eq. (2.68). 

When considering low temperature physics in a background magnetic held, -C 
|H|, it is natural to choose a renormalization point /i = 0(|i3|^/^), as this is the relevant 
scale which cuts oh long range huctuations in the charged SYM helds. We will employ 
two choices for the renormalization point. One choice is /i = 1/T, with L the AdS 
curvature scale (discussed below); this choice is computationally convenient but not 
physically signihcant. We will also report and discuss results with /i = For 

^°Note that metric variation of some renormalized EM action (whose separation 

from the total action would not be well-dehned). Rather, eq. (2.12) is simply defining T’g^(^) as the 
classical EM stress-energy tensor multiplied by the scale-dependent inverse EM coupling. 




later convenience, we define abbreviations for the (rescaled) energy density evaluated 
at these two renormalization points, 


£l = £(1/L), £8=£(|e|'''^). 


(2.17) 


2.2 Holographic description 

The holographic description of SYM states, within our sector of interest, in the limit 
of large Nc and large ‘t Hooft coupling A, is given by classical Einstein-Maxwell theory 
on 5-dimensional spacetimes which are asymptotically AdSs [42], The 5D bulk action 


IS 



(2.18) 


with G 5 = the 5D Newton gravitational constant, A = — 6 /L^ the cosmological 

constant, and L the AdS curvature scale.Setting to zero the variation of the action 
with respect to the metric gives the Einstein equation. 


Rkl + (A — ^R) Gkl — {FkmFl^ — \Gkl FmnF^^) , (2-19) 


while varying the bulk gauge held (with Fmn = '^m^n — '^n^m) gives the usual 
sourceless Maxwell equation, VrF^^ = 0. 


A 5D Chern-Simons term, A A F A F, could be added to the action (2.18) and 


would appear with a known coefficient in a consistent truncation of lOD supergravity. 
(See, for example, refs. [40, 42].) However, as stated above, in this paper we consider 
solutions with non-zero chemical potential /i or non-zero magnetic held B, but not both 
/i and B non-zero. For such solutions, the Chern-Simons term makes no contribution 
to the dynamics and hence may be neglected. 

As usual in holography, the expectation value (T"^(x)) of the stress-energy tensor 
is determined by the sub leading near-boundary behavior of the 5D metric Gmn- The 
leading near-boundary behavior of the bulk gauge held Am will be hxed by our chosen 
external U{1) gauge held (2.2), while the expectation value (j"(x)) of the f/(l) current 

^^The coefficient of the Maxwell action may, of course, be set to an arbitrary value by suitably 
rescaling the bulk gauge field Am- However, as the on-shell variation of the gravitational action with 
respect to the boundary value of the gauge field defines the associated current, such rescaling changes 
the normalization of the U{1) current in the holographic description. It will be seen below that 
the coefficient of the Maxwell term in our action (2.18) is correctly chosen so that the 17(1) current 
normalization is consistent with our previous charge assignments. If charge assignments are chosen, 
for example, to be larger by a factor of -s/S, then either the Maxwell term in the action (2.18) must 
be multiplied by a factor of 3, or else one must regard the boundary value of the bulk gauge field as 
equaling -s/S times the QFT gauge field (and the charge density in the bulk theory as equal to the 
QFT charge density divided by v^), as was done in ref. [40]. 
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density is determined by the subleading near-boundary behavior of the bulk gauge held. 
The precise relations will be shown below. 

Following ref. [18] , we choose a coordinate ansatz, based on generalized Eddington- 
Finklestein (EF) coordinates, which is natural for gravitational infall problems. The 
metric has the general form 

ds^ = — gajsix^r) dx°‘dx^ — 2wa{x) dx'^dr, (2.20) 

where r is the bulk radial coordinate and x = {a;"}, a = 0, ••■,3, denotes the four 
remaining spacetime coordinates. The spacetime boundary lies at r = oo; the {a:“} 
may be regarded as coordinates on the spacetime boundary where the dual held theory 
“lives”. Curves of varying r, with x held hxed, are radially infalling null geodesics, 
affinely parameterized by r. The one-form w = Wa dx°' (which is assumed to be time¬ 
like) depends only on x, not on r. These infalling coordinates remain regular across 
future null horizons. 

The form of the ansatz (2.20) remains invariant under r-independent diheomor- 
phisms, 

x“ ^ x“ = /"(x), (2.21) 

as well as radial shifts (with arbitrary x dependence), 

r ^ f = r + A(x). ( 2 . 22 ) 


We use the diheomorphism freedom (2.21) to transform the timelike one-form w to the 
standard form —dx^ (or Wa = —5°). Our procedure for dealing with the radial shift 
invariance ( 2 . 22 ) is discussed below in subsection 2 . 6 . 

We are interested in geometries which, at large r, asymptotically approach (the 
Poincare patch of) AdSs. This will be the case if gag^x^r) approaches as r — )■ cxd, 
with gaji = diag(—1,1,1,1) the usual Minkowski metric tensor. Demanding that the 
metric and bulk gauge held satisfy the Einstein-Maxwell equations, one may derive 
the near-boundary asymptotic behavior of the helds. Using radial gauge, Ar = 0, for 
the bulk gauge held, and a suitable choice of the radial shift ( 2 . 22 ) (which eliminates 
0 (l/r) terms in gap), one hnds that for solutions of interest, the metric and gauge held 
have asymptotic expansions of the form 


gapix,r) ~ gap + 



(LVr)^ + 0[(LVr)^], 


Aa{x,r) ~ A“*(x) -F A^^\x) {L^/rf + 0[{L^/r)^]. 


(2.23a) 

(2.23b) 
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The coefficient of the logarithmic term in the metric is only non-zero when there 
is an external EM held, 


C = F^urp 






(2.24) 


For a constant magnetic held in the direction, diag(-|-2, -f-l, -|-1, —2). 

The subleading asymptotic coefficients {x) and Aa '^ {x) cannot be determined solely 
from a near-boundary analysis of the held equations, and depend on the form of the so¬ 
lution throughout the bulk. However, asymptotic analysis does show that ~ 

— ^Fqj.Fq'^. The subleading metric coefficients g^pix) and encode the expec¬ 

tation value of the SYM stress-energy tensor [43, 44]. The appropriate holographic 
relation is 

= K I gl^^ - tr + [ln(/iT) + C] hjfj |, (2.25) 


where 


qA) — qA) 4- ^ n 
4 / 11 / — ' 4 


1 


gta + 4 ^00 


+ i > (2-26) 


K = L^/(47rG7v) = {N^ — l)/(27r^), and C is an arbitrary renormalization-scheme 
dependent constant.We adopt a specihc value. 


C = -|, (2.27) 

which will make the subsequent explicit expression (2.42a) for the energy density as 

simple as possible. 

Inserting expression (2.24) into relation (2.26) shows that the coefficient of 
the renormalization point dependent part of the holographic SYM stress-energy, is 
proportional to the classical EM stress-energy tensor, 

- I g^p F^.F ^'', (2.28) 

or = I diag(-|-l, -|-1, -|-1, —1) for a constant magnetic held in the x^ direc¬ 
tion. Using the above relations, one also hnds that tr(^*^"^^) = A . Since is 

traceless, the holographic relation (2.25) yields the stress-energy trace 

(T^) = -3«:tr(^(")) = -\nF^,F^\ (2.29) 

FefFerman-Graham (FG) coordinates, for which ds^ = {L?/p^) \gap{x, p) dx°‘ , one has 

dapix, p) ~ Pap + Idap (x) + In + 0(p® Inp) as p — ?> 0. Eq. (2.26) gives the relation between 
the subleading asymptotic metric coefficients in our infalling EF coordinates and FG coordinates. 

^^To perform the required holographic renormalization one must add a counterterm depending 
logarithmically on the UV cutoff. (See, for example, refs. [44-46].) As always, such a logarithmic 
counterterm comes with an inevitable finite ambiguity. 
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or {T"„) = in a constant magnetic field, in agreement with the earlier field 

theory result (2.5). Similarly, the renormalization point dependence of the stress-energy 
(2.25) coincides with the QFT result (2.14).^"^ 

Finally, the subleading asymptotic coefficient Aa\x) for the bulk gauge held en¬ 
codes the f/(l) current density. One hnds 

{j,) = 2kAI^K (2.30) 


2.3 Symmetry specialization 

As noted earlier, we are interested in studying solutions of Einstein-Maxwell theory 
which are spatially homogeneous. This implies that all metric functions depend only 
on and r. The arbitrary function A in the residual radial shift diffeomorphism (2.22) 
will depend only on Henceforth, for convenience, we will use n as a synonym for 

n is a null time coordinate. (In other words, v = const, surfaces are null slices of 
the geometry.) At the boundary, v coincides with the time t of the dual held theory. 

We also impose invariance under 0(2) rotations in the x^-x‘^ plane. This implies 
that only the goo, gos, gos, and gn = g 22 components of ga/s are non-zero. Our Einstein- 
Maxwell theory (without a Chern-Simons term) is also invariant under spatial parity, 
or x^ —>■ —x^ rehections, and for simplicity we will also impose parity invariance. This 
requires the vanishing of gos- 

For the bulk gauge held, the choice of radial gauge, A^ = 0, plus our imposed 
symmetries imply that 

A«(a:,r) = A“*(a:)-0(n,r)h°. (2.31) 

The corresponding bulk held strength, which is what appears in the held equations, 
can have a constant {x and r independent) magnetic held plus a radial electric held, 

Fi 2 {x,r) = B, For{x,r) = dr(j){v,r) =-£{v,r), (2.32) 


with all other components vanishing. 

As in ref. [18], it is convenient to rename the non-vanishing metric components as 

^f?oo = -2A, = = (2.33) 

where A, B, and S are functions of v and r. The resulting line element is 

ds^ = 2dv [dr — A{v, r) dv] + S(n, r)^ (dx^ -|- dy^) + . (2.34) 

in ref. [40], one can also use a comparison of holographic and QFT evaluations of the U{1) 
anomaly to confirm that the C/(l) current normalizations are consistent. 
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Henceforth, A will always denote the metric function multiplying dv^ (times —1/2), 
not the bulk gauge held. The function S is the spatial scale factor (with dx dy dz 
the spatial volume element), while B characterizes the spatial anisotropy (which should 
not be confused with the magnetic held amplitude B). 

The radial derivative dr is a directional derivative along infalling radial null geodesics. 
It proves convenient to dehne a corresponding directional derivative along outward ra¬ 
dial null geodesics, 

d+ = dr + A{v, r) dr . (2.35) 

The held equations which result from varying the action (2.18), inserting the above 
symmetry specializations, and re-expressing n-derivatives in terms of the (i+ modihed 
time derivative (2.35), take a remarkably compact form. The Einstein equations are: 

E" + |(H')^E = 0, (2.36a) 

- 6(E7s 2) d+S + I B'd+B = + I , (2.36b) 

(d+H)' + KS'/S) d+H + \B' (d+S)/S = -\B‘^L^ , (2.36c) 

(d+S)7S + 2(S7 s 7 d+S = -\B^L^ - \E‘^L^ + 2/^^ , (2.36d) 

d+(d+S) - A' (d+S) + lE (d+Bf = 0, (2.36e) 

where primes denote radial derivatives, h' = drh. As discussed in ref. [18], the 
anisotropy function B encodes the essential propagating degrees of freedom. The func¬ 
tions S and A may be regarded as auxiliary helds, determined by solving eqns. (2.36a) 
and (2.36b) using data on a single time slice. Information about the time evolution of 
B is contained in equation (2.36c). Equations (2.36d) and (2.36e) may be viewed as 
boundary value constraints — if they hold at one value of r, then the other equations 
ensure that these equations hold at all values of r. 

Maxwell’s equations reduce to the statements that neither the magnetic held H, 
nor the radial electric hux density £ have any radial or temporal variation. In other 
words, B = const., as already indicated in (2.32), and 

T(n,r) = pLS“^(n,r), (2.37) 

for some constant p which, from eqs. (2.30)~(2.32) plus (2.40) below, one sees is precisely 
the U{1) charge density (rescaled by k), 

{f)=Kp. (2.38) 

The form (2.37) of the radial electric held simply rehects Gauss’ law in 4-1-1 dimensions, 
combined with charge conservation and spatial translation invariance, which imply that 
p cannot have any temporal or spatial variation. 
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The bulk gauge field Am does not appear in the field equations (except via the field 
strength), but one may choose to regard Am as satisfying the radial gauge condition, 
Ar = 0, plus the condition that the time component A^ vanish at the horizon. This 
fixes the residual r-independent gauge freedom which remains after imposing radial 
gauge. With these choices, the chemical potential fi is the boundary value of A^ in 
the late time {v —)■ cxd) equilibrium limit. Equivalently (in radial gauge), the chemical 
potential /r equals the difference between the boundary and horizon values of Ay, in the 
equilibrium geometry. This coincides with the line integral of the radial electric field 
from horizon to boundary, 

poo 

= lim = / dr £{oo,r), (2.39) 

’’ Jri, 

which gives the work needed to move a unit charge from the boundary to the horizon. As 
usual, the charge density and the chemical potential are thermodynamically conjugate. 
One may consider the chemical potential /x to be a function of the (rescaled) charge 
density p, or vice-versa. The choice of perspective (“canonical” vs. “grand canonical”) 
has no bearing on the dynamics. 

2.4 Asymptotic analysis 

Asymptotic analysis of these equations is straightforward. We impose a fiat boundary 
geometry with the requirement that limr_>.oo 9 ap{x,r) = implying 

Wm. {L/rY A{v,r) = \ , lim (L/r) S(u, r) = 1, lim i?(u, r) = 0 , (2.40) 

r^oo ^ r—^oo r^oo 

for our renamed metric functions. Solutions to Einstein’s equations (2.36) with this 
leading behavior may be systematically expanded in integer powers of 1/r and (for 
non-zero magnetic field) logarithms of r. One finds: 

S(u,r) ~ L-^r + A(u)] + 0[{L/rY £], (2.41a) 

A{v, r) ~ A(u)]^ - dyX{v) 

+ [a4 - In f ] (L/r)^ 

- [204 A(u) + IB^ X{v) (1 - 2 In £)] {L/rY + 0[{L/rY In £], (2.41b) 

B{v,r) ~ L* [bYv) + iBAnj^] {L/rY 

+ [L^ dyhYv) - AbYv) A(u) + IB"^ A(u) (1 - 4 In f)] {L/rY 

+ 0[(L/r)^ln£], (2.41c) 

The constant 04 and the function b^v) cannot be determined just using asymptotic 
analysis, and the radial shift X{v) is completely arbitrary. The coefficient 04 encodes 
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the energy density which, due to homogeneity, cannot vary in time, while & 4 (n) encodes 
the anisotropy in the spatial stress. Using the holographic relation (2.25) and our 
convention (2.27) for dehning the stress-energy tensor, one hnds 

(T°°) = K (-§04 + I In fiL) , (2.42a) 

(T^i) = (T22) = (-la 4 + b^-lB^ + l B^ In/iL) , (2.42b) 

(T33) = K (-|a 4 - 264 - ^ B'^ In/xL) . (2.42c) 

2.5 Scaling relations 

Consider independent rescaling of the boundary and radial coordinates, 

X = ax, r = a!“^ 7 ^r, (2.43) 

with a and 7 arbitrary positive numbers. If the metric functions {A, T,,B} satisfy the 
Einstein equations (2.36), with asymptotic behavior (2.40), then the rescaled metric 
functions 

B{x, r) = B{x{x),r{j)) , (2.44a) 

S(5;, r) = (a/ 7 ) S(a;(5;), r(r)), (2.44b) 

A{x,r) ={a/'-^YA{x{x),r{j)), (2.44c) 

also satisfy the Einstein equations (and our asymptotic conditions) with rescaled pa¬ 
rameters 

L = L ^ B = B, p = a^ p. (2-45) 

The subleading asymptotic coefficients 04 and 64 become 

04 = a^ [04 — I ln( 7 /a)] , 64 = a'^ [ 64 -|-| ln( 7 /a)] . (2.46) 

Using the holographic relation (2.42) for the stress-energy expectation, one hnds that 

f^'^(^) = a^T'^^(/i), (2.47) 

with a rescaled renormalization point J1 = a p. 

If a = 7 , then these transformations are just a trivial rescaling of all quantities 
according to their dimension. But transformations with a 7 ^ 7 are non-trivial and scale 
bulk and boundary quantities by different amounts. In particular, transformations with 
a = 1 but 7 7 ^ 1 rescale the AdS curvature scale L without affecting the boundary 
coordinates or boundary parameters {B, p, or p), showing that the value of L has no 
physical signihcance (in the large Nc, large A limit for which classical gravity provides 
the dual description). This illustrates, explicitly, the independence of the boundary 
held theory on the AdS curvature scale L. 
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Figure 1. With a generic choice of the radial shift X{v) (left panel), the radial position of 
the horizon will change with time. It may be kept fixed (right panel) with a suitable choice 
of X{v). 

2.6 Apparent horizon 

With a non-zero homogeneous energy density, the dual geometries of interest will have 
an apparent horizon at some radial position, r = rh(n) [18]. Since we are investigating 
non-equilibrium dynamics, one might expect the horizon position to change significantly 
before ultimately settling down as the system equilibrates. However, as illustrated in 
£g. 1 it is possible, and very convenient, to use the residual radial shift diffeomorphism 
freedom (2.22) to place the apparent horizon at a hxed radial position, 


rh{v)=ri,. (2.48) 

A short exercise [18] shows that the condition for an apparent horizon to be present 
at r = fh is that this location be a zero of the modihed time derivative of the spatial 
scale factor, 

d+E|^^^_^ = 0. (2.49) 

This condition serves to £x the radial shift A(n). It is convenient to regard this condition 
as a combination of a constraint on initial data (implemented by finding the radial 
shift A(no) which is needed to satisfy (2.49) at some initial time Uq) together with the 
requirement that the horizon be time-independent, dr^^/dv = 0, which requires that 
the time derivative of d+S vanish at the apparent horizon. Evaluating this condition, 
and using the Einstein equations (2.36d, 2.36e) to simplify the result, determines the 
value of the metric function A at the horizon, 

■4 W =-1 . (2.50) 

For the metric to be non-singular on and outside the apparent horizon, the spatial scale 
factor S must be non-vanishing for r > fh. 
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2.7 Equilibrium solutions 

Given some initial non-equilibrium state of the system, the dynamical evolution should 
asymptotically approach a thermal equilibrium state. In the gravitational description, 
this implies that the geometry should, at late times, approach some static black brane 
solution. The specihc black brane solution will depend on the values of the conserved 
energy and charge densities in the chosen initial state, and on the value of the back¬ 
ground magnetic held. 


Schwarzschild 

For initial states with vanishing charge density and magnetic held, the bulk geometry 
will equilibrate to the 5D AdS-Schwarzschild black brane solution. A standard form of 
this metric is 


df2 ~2 


{i = 1, 2, 3), with 


ds =-U{f)dt + 


U(f) ^ f! _ !!A 

^ ’ L2 P 


i\2 


(2.51) 


(2.52) 


The radial coordinate r should not be confused with our Eddington-Finklestein coor¬ 
dinate r. The zero of U{r) determines the horizon location, 

Th = , (2.53) 

and the horizon temperature [given by (27r)“^ times the surface gravity at the horizon] 
is proportional to the horizon radius. 


TrTh = r^L ^ /L . 


(2.54) 


In our infalling EF coordinates (2.34), this AdS-Schwarzschild solution is described by 

S(r) = (r+A)/L, A{r) = , 5(r) = 0. (2.55) 

Using the holographic relation (2.42), one sees that the parameter m is related to the 
(rescaled) equilibrium energy density £ = /k via 

£=|mL-U (2.56) 

Reissner-Nordstrom 


If the initial state has a non-zero charge density but vanishing magnetic held, then the 
bulk geometry will equilibrate to a 5D Reissner-Nordstrom (RN) black brane [42]. This 
metric may be written in the form (2.51), with 

+ (2.57) 
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Figure 2. The one parameter family of non-extremal equilibrium Reissner-Nordstrom 
charged black brane solutions (solid line) in the plane of el\p\^l‘^ vs. {'kT)^/\ p\^/^. Also 
shown are the high and low temperature asymptotic forms (dashed lines). In the high tem¬ 
perature regime, ttT » the curve approaches the Schwarzschild result, e = |(7rT)^. 

In the low temperature (or near extremal) regime, vrT <C the charge density p ~ 

Pm„[l - and e/p-/’ ~ (1)-/’ + 1 (3)--/’ {er)V--/3. 


The charge density p of the Reissner-Nordstrom brane is bounded from above by the 
extremal charge density Pmax, given by 


{Pn 


rii\4 _ 4 3 

^±j j 2 * 


(2.58) 


The relation (2.56) between the energy density and the mass parameter m is unchanged. 
Hence, the extremal charge density pmax = | 

It is convenient to express p in terms of the fraction x of the extremal charge 
density, 

X = p j Pmax • (2.59) 

The horizon radius bh is given by the outermost positive root of f/(r); explicitly. 


rh/T = 


{—x^ + i\/l — x'^Y^^ + {—x^ — Wl — 


11/2 


(2.60) 


The horizon radius (divided by L) varies from down to (|m)^/'^ as x varies from 
0 to 1. The horizon temperature is given by 


vrThL = {fjL) [1 - (|m)3/2 . (2.61) 


The horizon temperature decreases with increasing charge density, and vanishes as the 
charge density approaches the extremal value (or a: —>■ 1 ). 
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From the perspective of the dual held theory, for any given value of the charge 
density there is a lower bound on the energy density, £min(p), which must be a mono- 
tonically increasing (and convex) function of p. This implies that for any given value of 
the energy density, there will be a maximum charge density, corresponding to a ground 
state with vanishing temperature. The equilibrium chemical potential /r, thermody¬ 
namically conjugate to the charge density p, is given by 

p=\p{L^lf^f. (2.62) 

Physically distinct non-extremal solutions may be labeled by the value of one dimen¬ 
sionless ratio such as e/\p\^f^ [or (7rT)^/|p|^/^ or e/(7rT)"^, etc.]. Figure 2 shows a log-log 
plot of the curve representing these solutions in the plane of e/\p\^/‘^ and 

In our infalling EF coordinates, the RN black-brane solution is described by 

S(r) = (r+A)/L, A{r) = - \m^{r)-^ + \p^L^T.{r)-\ (2.63) 

and B{r) = 0. 

Magnetic branes 

When the magnetic held is non-zero, the bulk geometry will equilibrate to a stationary 
magnetic black brane solution. These solutions are not known analytically, but have 
been studied numerically [40, 43]. In our infalling coordinates, the solutions satisfy the 
static specialization of eqs. (2.36).^^ The near-boundary behavior of asymptotically 
AdSs solutions is given by the expansions (2.41) (but with no time dependence). 

The extremal, zero-temperature magnetic brane solution interpolates smoothly be¬ 
tween AdSs X near the horizon and AdSs near the boundary. In our infalling 
coordinates, a series in fractional powers of 5r = r—rh describes deviations from the 


AdSs X 1^^ geometry near the horizon, 

A{r) = |((5r/L)^ [l + h + 0{rf , (2.65a) 

S(r) = {BL6ry/^ [l - p + 0{r]^ 6r^"^)] , (2.65b) 

B{r) = -|ln[275r7(i32L8)] - T( 13 + 27 ) + 0(7 , (2.65c) 

^®The resulting equations may be written explicitly as: 

(AE^)' E-3 = E-4 + + 4 / 1,2 ^ (2.64a) 

(TE'E^)' E-3 = J.-4 _ 1 ^ 2^2 ^j , (2.64b) 

[AB'Y?)' E-3 = E"^ . (2.64c) 
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with 7 = —1 + The constant r] cannot be determined from a near-horizon 

analysis and mnst be snitably adjnsted after integrating eqs. (2.64) to obtain the desired 
bonndary geometry. There is a single extremal magnetic brane solntion, modnlo the 
rescaling transformations (2.43)-(2.45) (which relate solntions with any non-zero valnes 
of the magnetic held B and cnrvatnre scale L), and radial shift diffeomorphisms (2.22). 

For non-extremal solntions (with non-zero B bnt vanishing p), metric fnnctions 
near the horizon have power series expansions in 6r = r—rh of the form 


A{r) = oo 5r L~^ [l — (1 — | B^L'^Sq^) -|- 0(5r^)] , (2.66a) 

S(^") = So/l3 [l-f (2 - I B'^L^Sq^) Oq ^ 6r + 0(5r^)] , (2.66b) 

B{r) = 2 In /d — I B‘^L^Sq‘^ 5r -|- 0{5r'^). (2.66c) 

The coefficient oq is proportional to the horizon temperatnre, 

T = ao/(27rL2). (2.67) 


The other two nndetermined constants, Sq and /d, which control the horizon valnes of 
the spatial scale factor and the anisotropy fnnction, mnst be snitably adjnsted after 
integrating eqs. (2.64) to select solntions which have the desired near-bonndary behavior 
(with an isotropic bonndary metric). If B^ ^ T'^, then the resnlting magnetic brane 
geometry is a small pertnrbation away from the Schwarzschild black brane (2.55), while 
if B^ ^ then the geometry may be regarded as interpolating between the BTZ black 
brane (x M^) near the horizon and AdSs near the bonndary [40]. 

There is a one parameter family of non-extremal solntions, modnlo the rescaling 
transformations (2.43)-(2.45) (and radial shift diffeomorphisms). Physically distinct 
solntions may be labeled by the valne of the dimensionless ratio [or {itTY/B'^ 

or £g/(7rT)‘^, etc.]. The left panel of hgnre 3 shows a log-log plot of onr nnmerically 
determined cnrve representing these solntions in the plane of Sb/B'^ and [ttTY/B'^. 
Extrapolating onr lowest temperatnre nnmerical resnlts to zero temperatnre, we hnd 
estimates of 

Co ^0.18, id*(p) ^ 2.0/i^ (2.68) 

for the coefficient cq or the eqnivalent hdncial scale B* dehned by eqs. (2.15) and (2.16). 

If one chooses to measnre energy density and magnetic held in nnits set by the 
cnrvatnre scale L, then one may traverse the one-parameter family of magnetic brane 
solntions by varying \B\L‘^ for a hxed valne of (or vice versa). The holographic 
relation (2.42) shows that these cnrvatnre scale dependent qnantities are related to the 
intrinsic dimensionless parameter EbIB"^ via 

i = + (2.69) 
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Figure 3. Left: The one parameter family of non-extremal equilibrium magnetic brane 
solutions in the plane of vs. /B'^. Also shown are the high and low temperature 

asymptotic forms (dashed lines). For high temperatures, vrT the curve approaches 

the Schwarzschild result e = |(7rr)^. For low temperatures, ttT <C the form ebIB'^ ~ 

Cl + C 2 {-kTY/\B\ provides a good fit to our data for ci = 0.35 and C 2 = 0.20. Right: The 
relation between the intrinsic parameter labeling magnetic brane solutions and the 

value of the magnetic field in curvature scale units, \B\L?‘, for two different fixed values of the 
curvature scale energy density, elL^ = ±0.75. 

This relation between and \B\L‘^, for two different fixed values of the curvature 

scale energy density = ±0.75, is plotted in the right panel of fig. 3. Note that two 
different values of \B\L‘^ yield the same value of Sb/B"^ (and hence describe the same 
physical solution) when eg > 0 . 

For these non-extremal magnetic brane solutions, the anisotropy function B{r) 
increases (and the scale factor S decreases) smoothly as one moves inward from the 
boundary toward the horizon. Figure 4 (left) plots the resulting anisotropy function 
B[r) for several values of the magnetic held when the energy density at the scale 1/L 
is held hxed, SlL^ = |. (From eq. (2.42a), this is the same as hxing the asymptotic 
coefficient 04 /^^ ~ horizon temperatures for this series of solutions, in order 

of increasing magnetic held, are given by ttTL = 0.873, 0.806, 0.879, 1.103, and 1.347. 
From the hgure, one may see that the horizon value of the anisotropy function is not a 
monotonic function of magnetic held (for hxed a^L"^). 

The right panel of hg. 4 shows a similar set of solutions with increasing magnetic 
held, but now with the energy density at the scale \B\^/'^ held hxed, Eb = 8 L“^. With 
the physical parameter Eb held hxed, the horizon value of the anisotropy function is 
now monotonically increasing with magnetic held. The temperatures of these solutions 
(in order of increasing B) are given by ttTL = 1.807, 1.797, 1.738, 1.620, and 1.433. 


21 















B 


B 




Figure 4. The anisotropy function B{r) as a function of inverse radius u = 1/r for equilibrium 
magnetic brane solutions with different values of the magnetic field. Left panel: Solutions 
at fixed energy density el = 0.75 with BL^ = 1.0, 1.5, 2.5, 3.5, and 4.5. Right panel: 
Solutions at fixed energy density eb = 8.0 with BL?' = 0.0, 1.0, 2.0, 3.0, and 4.0. In all 
cases, the radial shift A has been suitably adjusted to fix the horizon radius at u = 1. Note 
that the horizon value of the anisotropy function is not a monotonic function of magnetic 
field at fixed el, but is monotonic when eg is held fixed. 

3 Techniques 

3.1 Computational strategy 

We apply the computational strategy presented in ref. [18] to our case of homogeneous 
isotropization in Einstein-Maxwell theory. For convenience, we choose units in which 
the AdS curvature scale L = 1 . 

Required initial data, on some v = vo time slice, consists of an initial choice for 
the anisotropy function B{vo,r) and the radial shift A(uo )5 along with chosen (time 
independent) values of the energy density e, charge density p, and magnetic field B. 
As noted above, the holographic relation (2.42a) shows that fixing the energy density 
e: at the scale /i = 1/L is equivalent to fixing the asymptotic coefficient a^. Our choices 
for the initial anisotropy function will be detailed below in subsection 3.3. 

Given a set of initial data, the linear second order radial ordinary differential equa¬ 
tion (ODE) (2.36a) may be integrated to hnd the spatial scale factor S(uo,r). The 
two leading terms in the asymptotic behavior (2.41a) provide the integration constants 
needed to specify uniquely the desired solution. Next, one solves eq. (2.36d), a linear 
first order radial ODE for d+S. The near-boundary asymptotic behavior of this func¬ 
tion is d+S ~ |(r-|-A)^ -|- 04 -|- 0{r~^) and homogeneous solutions to eq. (2.36d) 

behave as near the boundary. Hence, the chosen value of the energy density e 
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uniquely specifies the desired solution. With B, E, and d+S determined on the Vq time 
slice, one next solves eq. (2.36c), a linear hrst order radial ODE for d+B. The desired 
asymptotic behavior of this function is d+B ~ — 2 h 4 r~^ + 0(r~^) while homogeneous 
solutions to eq. (2.36c) behave as near the boundary. So the needed integra¬ 

tion constant corresponds to requiring the absence of any such homogeneous solution. 
Finally, one solves the second order linear ODE (2.36b) to determine A{vQ,r). Homo¬ 
geneous solutions are linear or constant functions of r. From the asymptotic behavior 
(2.41b), one sees that the value of the radial shift A(uo) hxes the coefficient of the homo¬ 
geneous solution linear in r and provides one of the two needed boundary conditions. 
The second boundary condition, needed to fix the constant homogeneous solution, is 
provided by the horizon stationarity condition (2.50), which determines the value of A 
on the apparent horizon. 

Having solved for the modified time derivative d 4 .B{vo,r), and A{vo,r), one recon¬ 
structs the ordinary time derivative of the anisotropy function via 


drB{vo, r) = d+B{vo, r) - H(uo, r) drB{vQ, r). (3.1) 

The time derivative of the radial shift, 9t,A(uo)5 is extracted from the asymptotic behav¬ 
ior (2.41b) of A by evaluating the r —)■ cx) limit of H — |(r-|-A)^. These time derivatives 
of B and A provide the information needed to advance in time. Using a standard 
numerical integration scheme, one takes a small step forward in time, advancing v to 
Vo + Av for some timestep Av. Repeating this process, one progressively determines 
the metric functions on a sequence of equally spaced time slices, v = Vk = Vq + k Av. 
On each time slice, the asymptotic coefficient b^i^v), needed to determined the stress 
tensor (2.42), is extracted from the large r behavior of the anisotropy function B. (In 
the presence of a non-zero magnetic field, one extracts from the large r limit of 
a subtracted, rescaled version of the anisotropy function which removes the leading 
logarithmic piece in eq. (2.41c). This is detailed in the next subsection.) 

3.2 Numerical methods 

We use an inverted radial coordinate u = 1/r, and arbitrarily choose 

fh = l, (3.2) 

as our apparent horizon location. This makes our computational domain a fixed radial 
interval, 0 < u < Uh = 1- We use a 4th order Runge-Kutta method (described in 
ref. [18]) for time integration. This requires four integrations of our radial ODEs per 
time step, but yields much better accuracy, for a given timestep Av, than a lower order 
method. 
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To integrate the radial ODEs (2.36a-2.36d), we have used both traditional short- 
range hnite difference approximations, and spectral methods [47]. In the latter ap¬ 
proach, one implicitly represents the radial dependence of functions as a (truncated) 
series of Chebyshev polynomials. Explicitly, functions are represented by the list of 
their values on a discrete, finite collocation grid consisting of the points 



1 -|- cos 


TT k ■ 

M-l\ ’ 


/c = 0, • • ■ , M—1 , 


(3.3) 


and derivatives are represented by (dense) M x M matrices acting on the finite list of 
function values. The (truncated) spectral expansion converts each ODE into a straight¬ 
forward linear algebra problem. Boundary conditions are simply encoded into the first 
or last rows of the resulting matrix [47]. 

Although there are subtleties (described momentarily) in applying spectral methods 
to our problem, we have found the use of spectral methods to be clearly superior to 
finite difference approximations, yielding both faster computation and more accurate 
results. Using an M point discretization of the computation domain, short-range finite 
difference methods have errors which scale as an inverse power of M while spectral 
methods, in favorable cases, produce errors which fall exponentially with increasing M. 

Spectral methods presume that one is approximating functions which are regular 
and well-behaved on the computational domain. However, our radial ODEs have reg¬ 
ular singular points at m = 0 or r = cx) (due to the growth of the scale factor near 
the boundary), and our functions S, d+E, and A all diverge at the u = 0 endpoint. 
Therefore, we define subtracted functions in which the known singular near-boundary 
behavior is removed. To minimize loss of numerical precision in spectral approxima¬ 
tions of derivatives, it is also helpful to rescale the subtracted functions so that they 
do not vanish faster than linearly as u —)■ 0. If the magnetic field is zero, so no log¬ 
arithmic terms are present in the near-boundary behavior (2.41), then our subtracted 
functions are analytic in a neighborhood of the u G [0,1] radial interval and spectral 
methods converge exponentially. With a non-zero magnetic field, logarithmic terms 
are unavoidably present, showing that m = 0 is a branch point of the metric functions. 
This degrades convergence of the spectral series, leading to power-law convergence at a 
rate which depends on the behavior of the leading non-analyticity. Consequently, it is 
desirable to subtract logarithmic terms to as high an order as is practicable. We chose 
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to introduce subtracted/rescaled functions (denoted with a subscript ‘s’) via: 

E(r) = (r+A) + Ss(r), (3.4a) 

A{r) = |(r+A)^ — Inr — 2Xr~^ + 3A^r“"^ — 4A^r“^] + As{r). (3.4b) 

(i+S(r) = ^S(r)^ — Inr — 2Xr~^ + 3A^ — 4A^r“^] + ((i+S)s(r), (3.4c) 

d+B{r) = —Inr [r“^ — 3Ar“^ + 6A^r“^ — 10A^r“®] + r“^ {d+B)s{r). (3.4d) 

All our numerical work is performed using these subtracted/rescaled functions; we 
directly solve for S*, {d+J2)s, {d-^.B)s, and The expressions (3.4) are used to 

reconstruct the original functions when needed. We also use a subtracted/rescaled 
anisotropy function Bs{r), introduced via 

B{r) = Inr \r~^ — 4A r“® + lOA^r“®] + r~^Bs{r). (3.5) 

This removes leading logarithmic terms and introduces a convenient rescaling. Hence¬ 
forth, Bs will be referred to as the subtracted anisotropy function. 

In the above subtractions, the series in 1/r multiplying each logarithm are just 
truncated expansions of (r-|-A)“^ for A: = 2, 3 or 4. The choice to truncate these 
binomial series was arbitrary, but we found that our numerics were sufficiently well- 
behaved with the above subtractions. At higher orders in 1/r, additional terms appear 
which involve the asymptotic coefficient 64 and its (a-priori unknown) time derivatives, 
as well as higher powers of logarithms. 

3.3 Initial data 

As indicated above, one must select the value of the energy density (or asymptotic 
coefficient 04) and the initial value of the radial prohle of the anisotropy function, 
B{vo,r). And one must choose the value of the magnetic held B or charge density p. 
For the charged case, with vanishing magnetic held, physics can only depend on the 
dimensionless combination of charge and energy densities so a (positive) value 

of e may be chosen arbitrarily without loss of generality. Given a choice of e, possible 
values of the charge density p are limited by the extremality bound |p| < pmax = | 

For the initial anisotropy function, we chose to focus on Gaussian prohles. In the 
A = 0 frame, 

H(^;o,r) = (3.6) 

^®Note that (d+E)s ^ and likewise for these are just names for the sub¬ 

tracted/rescaled functions for d+E and d+B, respectively. 
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We investigate the dependence of results on the parameters of the Gaussian (amplitude 
A, width a, and mean position tq) in the hrst part of section Motivated by the fact 
that in our coordinates lines of varying r, at hxed v, are radially infalling geodesics, we 
will refer to this initial Gaussian as a “pulse” of initial anisotropy. 

For the magnetic case, as discussed above, the breaking of scale invariance implies 
the presence of logarithmic terms in the asymptotic behavior of the anisotropy function. 
We simply add the log terms shown in eq. (3.5) to the Gaussian (3.6). With vanishing 
radial shift, A = 0, our chosen initial anisotropy function takes the form 

B{vo, r)=A Inr . (3.7) 

Arguably, a more natural choice for the magnetic case initial data might be to add 
a Gaussian to the full equilibrium solution for the anisotropy function in the chosen 
magnetic held. This could be seen as nicely paralleling the charged case (in which 
the equilibrium solution has vanishing anisotropy). Nevertheless, we will stick with 
our somewhat arbitrary choice (3.7), which is an adjustable initial pulse added to 
the correct asymptotics. As will be seen, the Gaussian pulse will nearly always be 
the dominant portion of the deviation from equilibrium and the driving force of the 
resulting anisotropy in the boundary stress. We doubt that differing choices in the 
precise form of the slowly varying function to which the Gaussian pulse is added would 
impact, in any significant way, the characteristic equilibration times or other significant 
features of the results presented below. 

After choosing the initial anisotropy function (in the A = 0 frame) the initial value 
of the radial shift, A(no), is adjusted to £x the location of the apparent horizon, as 
discussed in section 2.6. 

It should be noted that, in all cases (charged, uncharged, magnetized) it is quite 
possible to select physically inconsistent initial data. This happens when the initial 
anisotropy, for a given energy density, is so large that no apparent horizon shields a 
coordinate singularity produced by a vanishing scale factor S. This sets a natural limit 
on the amplitude of the initial pulse which is meaningful to study. 

4 Results 

4.1 Neutral plasma 

Before presenting results for equilibration in charged or magnetized plasmas, we first 
discuss general features of the time evolution in the uncharged, unmagnetized case and 

i7For results from an exploration of a broader range of initial anisotropy profiles, in the case of 
vanishing charge density and magnetic field, see ref. [9]. 
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examine the sensitivity of results to specific features in the initial data. As noted above, 
we choose a Gaussian prohle (3.6) for the initial anisotropy function, with an adjustable 
amplitude, width, and mean position. Typical evolution of our subtracted/rescaled 
anisotropy function, Bs{v,u) = is shown on the left in hgure 5. One sees 

the initial pulse prohle on the back side of the plot at u = 0. The hgure clearly shows 
the inhuence of the pulse propagating outward and rehecting oh the boundary at u = 0 . 
The outgoing portion of the pulse essentially propagates along an outgoing radial null 
geodesic which, in our coordinates, near the boundary are 45° lines at constant values 
of M + u. The inhuence of the anisotropy pulse, after the rehection, largely falls through 
the horizon along an ingoing radial null geodesic which is instantaneous in our null time 
coordinate v. The asymptotic coefficient 64 , which equals the slope of at u = 0, 
controls the anisotropy in the stress tensor, 

AV = I (T^i) + I (T^^) - (T33) , (4.1) 

with AV/n = 3 64. Hence, the rehection of the pulse oh the boundary directly produces 
the pressure anisotropy AV in the boundary theory. The time dependence of the 
relative pressure anisotropy,^® dehned as AV/{Ke), is plotted on the right in hgure 5. 

As shown in the hgure, at late times the anisotropy function approaches zero, as 
required for equilibration to the isotropic Schwarzschild black brane solution. At suffi¬ 
ciently late times, when the departure from equilibrium is small, the evolution should 
be well described by a linearized approximation to the full nonlinear dynamics. The 
linearized dynamics of inhnitesimal perturbations away from equilibrium may be repre¬ 
sented as a sum of quasinormal modes (QNM), which are eigenfunctions of the linearized 
dynamics with complex frequencies, (j){t) = Re(Ae“*‘^*) with Imca < 0. The lowest 
quasinormal mode (for which —Imca is minimal) dominates the late time approach to 
equilibrium. For the Schwarzschild black brane, quasinormal mode frequencies have 
been previously evaluated by Starinets [48]. From the late time behavior of our full 
nonlinear evolution, it is straightforward to extract an estimate of the lowest quasinor¬ 
mal mode frequency. Comparing with the independent results of ref. [48] provides a 
useful test of the accuracy of our numerics. Fitting the late time (4 < < 25) 

portion of our calculated pressure anisotropy to a decaying, oscillating exponential, 
cos[(Rea;)u - 1 - 0], yields an estimate of the lowest QNM frequency u which 
agrees with ref. [48] to hve digits, u/{'kT) ^ 3.11946 — 2.74663 h 

We will see the same vanishing of the anisotropy function at late times for the 
case of charged plasmas, whose gravitational duals equilibrate to an isotropic Reissner- 
Nordstrom black brane solution. For the magnetic case, however, at late times there 

^®Note that, for unmagnetized plasma, the energy density is three times the average pressure, 
e = {TV = W. 
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Figure 5. Left: Rescaled anisotropy function Bg = u~^B, for a typical case of equilibration 
to the Schwarzschild black brane, as a function of inverse bulk radius u and time v. Initial 
pulse parameters are ^ = 5 x 10“^, ro = 4 and a = ^, with e = The (rescaled) energy 

density e is used to the set the scale for time. One sees that the effect of the initial Gaussian 
pulse propagates outward, essentially along an outgoing radial null geodesic, reflects off the 
boundary, and then largely falls through the horizon (along an ingoing radial null geodesic 
which is instantaneous in v). After one “bounce”, the anisotropy rapidly approaches zero. 
Right: The corresponding relative pressure anisotropy, /S.Vjne = — 2T^^)/«:e, 

induced in the boundary field theory, as a function of time. Note how the peaks of the pressure 
anisotropy correspond directly to the reflection of the anisotropy pulse off the boundary. 


is a non-zero profile for the anisotropy function, reflecting the spatial anisotropy of 
equilibrium magnetic brane solutions. 

We now turn to an examination of the dependence of the pressure anisotropy on the 
parameters of the initial Gaussian (3.6). Of particular interest will be the dependence 
of the response on the amplitude and position of the initial pulse. Less interesting is 
the dependence on the width of the pulse, which affects the duration of the reflection off 
the boundary (and also produces changes more naturally associated with the position 
of the pulse). 

In hgure 6, we compare the pressure anisotropies created by the pulse shown in 
£g. 5, and an otherwise identical pulse with half the amplitude. Visually, one sees that 
the smaller amplitude pulse produces roughly half the pressure anisotropy as does the 
larger pulse, but with a virtually identical time course. The peak pressure anisotropy 
(divided by energy density), for both pulses is over 4, signihcantly larger than unity. 
Hence both initial pulses represent large departures from equilibrium. Given the highly 
nonlinear nature of the Einstein equations, one might have expected to see clear signs of 
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Figure 6. Left: The pressure anisotropy bPj ke produced by the same initial pnlse shown 
in fig. 5 (blue curve), overlaid with the pressure anisotropy produced by an initial pulse with 
half the amplitude (purple curve). Halving the initial amplitude roughly halves the induced 
pressure anisotropy. Right; The “nonlinearity” (NL), defined as the difference between the 
pressure anisotropy produced by the larger pulse and twice the pressnre anisotropy produced 
by the half amplitude pulse. 

nonlinearity in the dependence of the pressnre anisotropy on the initial pulses. However, 
even for these pulses producing large departures from equilibrium, the amplitude of the 
peaks in the induced pressure anisotropy are nearly linear in the amplitude of the initial 
Gaussian pulse. 

The right hand panel of hgure 6 makes this comparison quantitative. This shows 
the nonlinearity (NL) dehned as the difference between the pressure anisotropy bV/ne 
of the larger initial pulse and twice the pressure anisotropy produced by the halved 
initial pulse. Compared to the pressure anisotropies themselves, the relative size of the 
nonlinearity is roughly one part in 10^. This suggests that the dynamics, as probed by 
these initial pulses, are surprisingly close to a linear dynamical system. 

In asymptotically AdS gravitational solutions, deviations of the geometry from 
that of pure AdS space necessarily vanish as one approaches the boundary. Hence, one 
might expect nearly linear dynamics to be evident for initial pulses which are localized 
sufficiently close to the boundary, while anticipating much larger nonlinearities for 
initial pulses localized closer to the horizon. To test this expectation, we used a very 
large Gaussian — the same initial Gaussian prohle which generated £g. 5 but with the 
amplitude increased^® by a factor of 40 — and then examined the resulting evolution 
when the mean position of the initial Gaussian was progressively shifted deeper into 
the bulk. Figure 7 compares the evolution for mean positions tq = {4,2,1, 4, 4 the 

^®For the chosen values of position and width of the Gaussian, plus energy density e = this 

amplitude is close to the upper limit set by demanding the existence of an apparent horizon. 
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Figure 7. Comparisons of initial anisotropy functions (left column), induced pressure 
anisotropy (middle column), and nonlinearity (right column) defined as in fig. 6, for a series 
of five Gaussian initial anisotropy functions differing only in their depth in the bulk. From 
top to bottom, the mean position of the initial pulse, in the A = 0 frame, is tq = {4, 2,1, 

In all cases the energy density is e = The plots in the hrst row come from the same 

initial data as in fig. 5, but with the amplitude increased by a factor of 40 {A = 0.02, tq = 4, 
a = ^). The left column shows the initial anisotropy function scaled by u~^ and plotted 
as a function of the inverse radial coordinate u, after adjusting the radial shift A to fix the 
apparent horizon at u = 1. 
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frame with radial shift A = 0. In all cases, the energy density was held hxed at a value 
(e = which puts the equilibrium horizon position at r = 1. In other words, the 

only change in the hve cases shown in £g. 7 is the radial position of the initial Gaussian 
anisotropy function (3.6), viewed as a function of r. 

Each row of fig. 7 displays results for one of these five cases. In each row, the 
left hand panel shows the initial anisotropy function, but plotted as a function of 
the inverse radial coordinate m, after adjusting the radial shift A to £x the apparent 
horizon at m = 1, as discussed in sec. 2.6. The middle panels show the resulting pressure 
anisotropy as a function of time, and the rightmost panels display the nonlinearity (NL), 
again dehned as the difference between the pressure anisotropy of the given initial pulse 
and twice the pressure anisotropy after halving the initial amplitude. 

From the middle column of plots, one sees that the magnitude of the pressure 
anisotropy decreases significantly as the initial pulse is moved deeper into the bulk. 
Moreover, both the time it takes for the effect of the pulse to reach the boundary, and 
the width of the resulting peaks in the pressure anisotropy, grow with increasing depth 
of the initial pulse. This reflects the usual holographic mapping between bulk and 
boundary: phenomena deeper in the bulk correspond to lower energy or longer time 
scales in the boundary held theory. 

Turning to the nonlinearity plots in the right hand column, one sees that the 
magnitude of the nonlinearity also decreases as the pulse moves deeper into the bulk. 
Dividing the peak nonlinearity by the peak pressure anisotropy gives a relative measure 
of nonlinearity. This is about 1 x 10“^ for the top row, 3 x 10“® for the middle row, 
and 3 x 10“® for the bottom row. So in this comparison, as the initial pulse is pushed 
deeper into the bulk, the relative nonlinearity decreases systematically. 

This comparison does not, however, imply that nonlinearities are never signihcant. 
The amplitude of the Gaussian pulse in the initial anisotropy function was kept hxed 
in hg. 7, resulting in a decreasing size of the induced pressure anisotropy as the pulse 
moves deeper into the bulk. While the hrst two rows of the hgure show pressure 
anisotropies which are large departures from equilibrium, AV/ne 3> 1, the hnal rows 
with AV/ne ^ 1 represent small departures from equilibrium. It is natural to ask 
what happens if one instead increases the amplitude as the pulse is moved into the 
bulk, so as to keep hxed the size of the induced pressure anisotropy. 

Such a comparison is shown in hg. 8, The upper row of the hgure shows the time 
dependence of the pressure anisotropy and the nonlinearity for the same pulse which 
generated hgs. 5 and 6 (M = 5 x ro = 4, a = |), while the lower row shows the 
pressure anisotropy and nonlinearity of a pulse with the same shape as the last row 
of hg. 7, but with larger amplitude {A = 2.5, ro = |, a = 4). The energy density 
remains hxed, e = \L~^. The peak pressure anisotropy is similar in the two cases. 
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Figure 8 . Top row: pressure anisotropy (left) and nonlinearity (right), defined as in fig. 6, 
as a function of time for the initial the pulse which created fig. 5. (^ = 5 x 10“^, rg = 4, 
a = ^). Bottom row: corresponding plots of pressure anisotropy and nonlinearity for a pulse 
located deeper in the bulk {A = 2.5, tq = j, cr = ^) with amplitude adjusted to produce a 
similar peak pressure anisotropy. In both cases the energy density is e = Substantial 

nonlinearity is present for this case, where an initial pulse deep in the bulk has sufficient 
amplitude to produce a large departure from equilibrium. 

and corresponds to a large departnre from eqnilibrinm. For the latter case of a pnlse 
deep in the bnlk, large enongh to indnce a far from eqnilibrinm pressnre anisotropy, the 
nonlinearity is signihcant, mnch larger than the previous examples. However, even for 
this case, the size of the nonlinearity relative to the peak pressure anisotropy is only 
about 5%, NL/{AV/Ke) ^ 0.05.^° 

The data shown in £g. 7 inspire several further questions. Looking down the middle 
column of the hgure, one sees that the onset of the response (i.e., the time of the hrst 
peak in the pressure anisotropy) increases as the initial pulse moves deeper into the 
bulk but seems, perhaps, to be approaching a maximum value. Is this really true, or 
can one craft initial data for which the onset of the response is much greater? As shown 

have also examined the level of nonlinearity by comparing the pressure anisotropy resulting 
from a sum of two different Gaussians to the sum of anisotropies induced by the individual Gaussians. 
The results were comparable to those discussed above and do not warrant separate discussion. 
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Figure 9. Initial anisotropy function (left), induced pressure anisotropy (middle), and non¬ 
linearity (right) for a narrow “deep pulse” ("4=i^)’'o = n)0' = ^) localized closer to the 
horizon than the deepest pulses of fig. 7. The energy density remains fixed at the same value, 
e = \L~^. Relative to the previous case of fig. 8, the induced pressure response has a delayed 
onset, but is otherwise very similar. 


on the left panels of the lower rows of the figure, when the average position ro of the 
Gaussian pulse is moved into the bulk, an increasingly large portion of the Gaussian 
ends up lying behind the apparent horizon. And when plotted as a function of our 
computational coordinate u = (f—A)“^ (with f the A = 0 frame radial coordinate), 
initial pulses with small values of tq are only moderately localized near the horizon — 
even though these pulses had constant widths when viewed as functions of r. As may 
be seen by comparing the left and middle columns of £g. 7, it is the leading edge of 
the anisotropy pulse (the region of near-maximal slope) which produces the hrst large 
response in the boundary anisotropy. Looking at the last two rows of the figure, one 
may question whether we are doing an adequate job exploring the response from initial 
disturbances which are localized close to the horizon. Will initial pulses which are more 
strongly localized near the horizon show significantly greater nonlinearity? 

Figures 9 and 10 show results of an effort to explore these questions. Fig. 9 shows 
the initial anisotropy function, along with the resulting pressure anisotropy and non¬ 
linearity, for a signihcantly narrower “deep pulse” {A = ’"o = if;Cr = ^)- A^^cl 

fig. 10 shows a 3D plot of the time dependent anisotropy function, plus the induced 
pressure anisotropy, for an extremely narrow deep pulse {A = ^, tq = 1, a = ^)-^^ 
The energy density in both cases remains fixed, £ = For the narrowest pulse, the 

amplitude ^ is near the upper limit which can be studied without destabilizing 
the horizon. In both figures 9 and 10, the peak pressure anisotropy AV/ne is large 
compared to unity, showing that these pulses are producing far from equilibrium initial 
states. 

As seen in these hgures, pulses which are more sharply localized very near the 
horizon do lead to a delayed onset in the resulting pressure anisotropy, occurring at 

^^This width is at the limit of what our numerics could do using a 240 point spectral grid. 















Figure 10. Rescaled anisotropy function Bg (left) as a function of inverted radius u and time 
u, for an extremely narrow pulse sitting at the horizon {A = ^, ro = 1, a = ^). Resulting 
pressure anisotropy (middle) and nonlinearity (right) as a function of time. Energy density 
e = The outward movement of the pulse toward the boundary clearly resembles the 

behavior of outgoing null geodesics originating very close to an event horizon, which can 
“hug” the horizon for extended periods before eventually escaping. 

~ 1.75 for the case of fig. 9 and ~ 3 for our narrowest pulse in fig. 10. But, 
within the range of pulse widths we have studied, the onset of the pressure anisotropy 
response is only delayed by a factor of 2-3 compared to the case of £g. 8. From the left 
panel of £g. 10, showing the time dependence of the (rescaled) anisotropy function 
one sees that the outward movement of the pulse toward the boundary resembles the 
behavior of outgoing null geodesics originating very close to an event horizon, which 
“hug” the horizon for extended periods before eventually escaping. One may wonder 
if the onset of the response in the pressure anisotropy could be delayed indehnitely by 
going to narrower and narrower initial pulses localized at the apparent horizon. We do 
not have a hrm analytic argument, but doubt that this is possible if one simultaneously 
demands that the amplitude of the response remain bounded away from zero.^^ 

The relative nonlinearity for the case of a narrow deep pulse shown in £g. 9 is quite 
small, about half a percent. But for the extremely narrow pulse with near maximal 
amplitude of £g. 10, the nonlinearity, relative to the pressure anisotropy, reaches the 
10% level. Linearization of the dynamics about equilibrium must provide an accurate 
approximation to the full nonlinear dynamics when the deviation of the geometry from 
the equilibrium Schwarzschild black brane solution is sufficiently small, as will be true 
at sufficiently late times. For asymptotically anti-de Sitter geometries, where metric 
functions have the asymptotic forms (2.41), this will also be the case for initial data 
involving perturbations localized sufficiently close to the boundary. (See refs. [8, 9, 49] 

^^This expectation reflects the diverging redshift of late emerging outgoing geodesics in the geomet¬ 
ric optics picture, and is consistent with the gapped spectrum of quasinormal mode frequencies for 
translationally invariant perturbations. 
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for related discussion.) It should be noted that reasonably good agreement between 
linearized dynamics and full nonlinear evolution was previously reported in ref. [8] and 
further explored in ref. [9] . In these works, the authors found agreement at a 20% level 
between the linearized and full dynamics for a variety of initial anisotropy prohles. Our 
results examining, more systematically, the dependence of the relative nonlinearity on 
the parameters of our initial Gaussian anisotropy function complement and extend this 
earlier work. Overall, despite prior knowledge of refs. [8, 9], we are still surprised by the 
small, often extremely small, levels of nonlinearity which we hnd even at early times 
when the induced pressure anisotropy is maximal, for initial perturbations localized 
deep in the bulk and producing large departures from equilibrium. 

4.2 Charged plasma 

We now turn to the equilibration of charged plasmas (by which we mean SYM plasmas 
with a non-zero density of the global U{1) conserved charge — not a plasma in which 
electromagnetic forces are included in the dynamics and Coulomb repulsion plays a 
signihcant role). As noted earlier in section 2.7, the bulk geometry should equilibrate to 
a non-singular Reissner-Nordstrom black brane solution provided the charge and energy 
densities satisfy the extremality bound (2.58), p < pmax = | However, for values of 
the charge density near Pmax, we hnd that sufficiently large initial metric perturbations 
can destroy the apparent horizon. We expect that such initial data are unphysical, not 
representing SYM initial states which could be produced by an operational procedure 
such as turning on time dependent external helds (which correspond to time dependent 
boundary data in the holographic description). In any case, we limit our attention to 
initial perturbations for which an apparent horizon is present at all times. We hnd that 
if one suitably decreases the amplitude of the initial departure from equilibrium while 
increasing the charge density, one can approach pmax from below while maintaining the 
existence of an apparent horizon. 

Figure 11 (left) compares the time dependence of the pressure anisotropy which 
results from initial data consisting of precisely the same Gaussian initial anisotropy 
function B{vQ,r) (in the A = 0 frame) and energy density as in hg. 5, and a charge 
density p equal to 0, 20, 40, 60, or 80% of the extremal density pmax- The immediately 
obvious qualitative result is that the hve diherent curves are so close together than they 
cannot be visually distinguished! Varying the charge density (at hxed energy density 
and hxed initial anisotropy function) has stunningly little impact on the subsequent 
time evolution. This is quantihed in the right panel of hg. 11 which plots the diher- 
ence in the pressure anisotropy /S.V/ne between the cases of p = 0.8pmax and p = 0. 
Gomparing the scales on the right and left hand plots, one sees that for this initial 
anisotropy function the sensitivity to the charge density is less than one part in 10^. 
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Figure 11. Left: time dependence of the pressure anisotropy (relative to ke ) for values of the 
charge density p which are 0, 20, 40, 60, or 80% of the extremal density pmax- The different 
curves are virtually indistinguishable. The initial anisotropy function B(vo,r) and energy 
density e are the same as in fig. 5. Right; difference in AV/ke between p = O.Spmax and 

p = 0. 

In the plots of fig. 11, we used the fourth root of the (rescaled) energy density, 
to set the scale for time. Since the energy density was held constant in the comparisons 
of £g. 11, this was a simple and convenient choice. For the hve cases shown in the 
hgure, e/{'KTY = 0.75, 0.76,0.79, 0.86, and 1.03 for p = 0, 20, 40, 60 and 80% of pmax, 
respectively. And, for comparison, the values of the equilibrium chemical potentials 
corresponding to these charge densities are given by p,/T = 0, 0.34, 0.73, 1.26, and 
2 .21, respectively. 

If the initial anisotropy pulse begins deeper in the bulk, then the sensitivity to 
the charge density is larger. Figure 12 shows a comparison of pressure anisotropies for 
differing charge densities, now using the deep pulse initial anisotropy function whose 
radial prohle has the shape shown in £g. 9.^^ As seen on the left panel, the amplitude of 
the response increases signihcantly as the charge density varies from 0 to 80% of Pmax- 
However, the time course of the equilibration (e.g., the times of the hrst or second peaks 
in the response, or the zero-crossing between these peaks) is only modestly affected, 
with changes of 3% percent or less. 

The fact that the sensitivity to the charge density is greatest for pulses close to 
the horizon is to be expected. In the equilibrium geometry (2.57), one sees that as 
the radius r increases from the horizon, the influence of the charge density decreases 

^^One might guess that the pulse used in the bottom row of hg. 8 would exhibit greater sensitivity 
to charge since it had a larger nonlinearity than the deep pulse of fig. 9. However, the latter pulse 
exhibits much greater sensitivity to charge. 
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Figure 12. Left; comparison of pressure anisotropies produced by different charge densities, 
up to 80% of extremality, when the initial anisotropy function is a “deep pulse” 
ro = cj = ^) with e = I Although the amplitude of the response grows, as shown, 
with increasing charge density, it is striking how little the time course of the response varies. 
Right: difference in pressure anisotropy /SV/ke between p = 0.8 pmax and p = 0. 

rapidly relative to the other terms in the metric. Only near the horizon, and close to 
extremality, does the charge density produce an 0(1) effect on the equilibrium geometry. 

At the beginning of this work, we expected that one interesting outcome would 
be information on the change in equilibration time produced by varying the plasma 
charge density. By “equilibration time”, we mean some rough but useful measure of 
when the departure from equilibrium is no longer substantial. To make this a bit more 
quantitative we adopt, somewhat arbitrarily, the criterion 

AP(t)/K£ < 0.1, (4.2) 

for all times t > teq, as indicating that the system is near equilibrium at time teq- 

Looking at the left panels of hgures 11 and 12, it is clear that the effect of the 
charge density on any reasonable measure of equilibration time can be summarized 
easily: there is very little effect! Even in Eg. 12, where the sensitivity to charge density 
is the largest we have found with horizon-preserving initial data,^'^ the time tpeak of 
the initial response peak and the approximate equilibration time teq both vary by less 
than 5%. 

Since figures 11 and 12 plot time in units set by the energy density, the high degree 
of insensitivity of the relaxation time course to the charge density seen in these hgures 

Achieving good numerical accuracy becomes increasingly difficult as one pushes toward extremal¬ 
ity, where the equilibrium solution bifurcates. Investigating very near extremal behavior more carefully 
is an interesting topic we leave to future work. 
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Figure 13. Comparison of the pressure anisotropies produced by two different charge densi¬ 
ties, 0% and 80% of extremality, when holding fixed the equilibrium temperature, ttF = 1 /L. 
The energy densities are 0.75 and 2.68, respectively. The initial anisotropy function is same as 
in fig. 5 (^ = 5 X 10“'^, ro = 4, cr = With time is plotted in units of inverse temperature, 
the relaxation time course shows negligible sensitivity to the charge density (although the 
amplitude of the response varies significantly). 

might lead one to think that the total energy density is playing a special role in setting 
the time scale of relaxation. Bnt it shonld be borne in mind that these hgures show 
comparisons in which, by design, both the initial anisotropy fnnction (in the A = 0 
frame) and the total energy density have been held fixed. Becanse the ratio of energy 
density to temperatnre (to the fonrth power) varies significantly with increasing charge 
density, it is natural to ask whether the degree of (in)sensitivity of the relaxation dy¬ 
namics to the charge density is substantially different if one holds fixed the equilibrium 
temperature instead of the energy density. Figure 13 shows such a comparison. Plotted 
are the pressure anisotropies resulting from the same initial anisotropy function used 
in hgs. 5 and 11, and charge densities of either 0 or 80% of extremality, but now with 
the energy density in either case suitably adjusted to £x the equilibrium temperature, 
ttT = 1/L. One again sees a significant change in the amplitude of the response with 
increasing charge density, but now with the temperature held hxed, increasing charge 
density decreases the amplitude of the response. Nevertheless, with time now plotted 
in units of (ttT)”^, one again sees negligible (~ 0 . 01 %) change in the time course of 
the equilibration dynamics as the charge density varies from zero and 80% of Pmax- 
Performing the same constant temperature response comparison using the deep 
pulse initial anisotropy function (whose radial prohle is shown in £g. 9), we hnd a 
larger — but still quite small — variation in the time course, approximately 2 %, as the 
charge density varies from zero to 80% of pmax- 

We have also examined the degree of nonlinearity in the above examples of equili- 
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Charged (p 7 ^ 0, 

6 = 0 ) 


P/Pmax 

IvaX/e^/^ 

A/ttT 

Linearized A/e^/'^ 

0.0 

3.35208 ±0.00004 -2.95144 ± 0.00013 

3.11946- 2.74663 i 

3.35207- 2.95150 i 

0.1 

3.34564 ± 0.00016 -2.95468 ± 0.00019 

3.11948 - 2.75763 i 

3.34568 - 2.95460 i 

0.2 

3.32624 ±0.00020 -2.96429 ± 0.00019 

3.13222 - 2.79139 i 

3.32630- 2.96444 i 

0.3 

3.29319 ± 0.00028 -2.98266 ± 0.00036 

3.14987- 2.85285 i 

3.29327- 2.98287 i 

0.4 

3.24572 ±0.00007 -3.01376 ± 0.00008 

3.17857- 2.9514H 

3.24574-3.01377 i 

0.5 

3.18366 ±0.00016 -3.06498 ± 0.00007 

3.22529 - 3.10506 i 

3.18370 - 3.0649H 

0.6 

3.11311 ± 0.00002 -3.15177 ± 0.00002 

3.31032-3.35142 i 

3.11311-3.15176i 

0.7 

3.07021 ± 0.00006 -3.29402 ± 0.00006 

3.50617- 3.76176 i 

3.07022 - 3.29399 i 

0.8 

3.11863 ± 0.00265 -3.41376 ± 0.00295 

3.99199 - 4.36977 i 

3.11848-3.42004 i 


Table 1. Lowest quasinormal mode frequency for charge densities ranging from 0 up to 80% 
of extremality. The second and third columns show results (with estimated uncertainties) 
for the real and imaginary part of the QNM frequency in units of obtained from fitting 
the late time behavior of our full nonlinear evolution. The fourth column shows these same 
results converted to units of ttT. (Fractional uncertainties are the same as in the preceding 
columns.) The rightmost column shows results from an independent analysis of the linearized 
small fluctuation equations by Janiszewski and Kaminski [50]. 

brating charged plasmas. The results for the relative size of the nonlinearity are quite 
similar to our earlier results for equilibrating uncharged plasmas. Because of this, we 
will refrain from presenting explicit nonlinearity plots for charged plasmas. 

Finally, as noted earlier, at sufficiently late times the relaxation must be accurately 
described by a superposition of quasinormal modes (eigenfunctions of the linearized 
dynamics about the equilibrium solution). Extracting the leading quasinormal mode 
frequency by htting the late time (4 < < 20) behavior of our calculated pressure 

anisotropy to a decaying, oscillating exponential, as described in the previous section, 
yields the results shown in table 1. The second and third columns (with uncertainties) 
show our estimates for the real and imaginary parts of the leading QNM frequency in 
units of while the fourth column (without uncertainties) shows our estimates con¬ 
verted to units of ttT. The rightmost column shows independent results of Janiszewski 
and Kaminski [50] obtained by analyzing the linearized small fluctuation equations 
about the RN black brane solution. The agreement is a satisfying conhrmation of our 
numerical accuracy. Interestingly, the imaginary part of the lowest QNM frequency 
varies by over 15% between p = 0 and p = O.Spmax- This is enormously larger than 
the part in 10^ sensitivity seen in £g. 11, and substantially bigger than the largest 
(~ 5%) sensitivity we found in the evolution time course with deep pulses, £g. 12. The 
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implications of these very differing sensitivities will be discussed further in section 5. 

4.3 Magnetized plasma 

We now present results of an analogous investigation of equilibration in plasmas (with 
vanishing charge density) in a homogeneous magnetic field B. Our discussion will 
parallel, as much as possible, the previous treatment of charged plasmas. But the 
breaking of scale invariance by the magnetic field, discussed in section 2 . 1 , produces 
some notable differences. One change, seen in section 2.4, is that the anisotropy function 
B{v, r) must now contain logarithmic terms in its near boundary behavior. We choose 
our initial anisotropy function to have the form (3.7) in which an adjustable Gaussian 
is added to the required leading logarithmic term. As in the previous discussion of 
charged plasmas, we will keep hxed the parameters of the Gaussian part of the initial 
anisotropy function B{vo, r) as we dial up the external magnetic field. We will also hold 
hxed the energy density dehned at a renormalization point fi = 1/L. As shown by the 
holographic relation (2.42), this is the same as holding hxed the asymptotic coefficient 
04 . In the following plots, axis labels involving energy density e, without any explicit 
indication of scale, will denote the energy density evaluated at the curvature scale, 
s{l/L) = El- Similarly, the pressure anisotropy AP should be understood as AP(1/L) 
unless otherwise indicated explicitly. 

Instead of keeping El hxed as the magnetic held B is varied, one could choose 
to hold hxed the energy density dehned at the scale set by the magnetic held, 
/i = Since the AdS curvature radius L is not a physical scale present in the dual 

QFT, hxing eb instead oi El is arguably more natural. However, for computational 
reasons it is easier to hold hxed El a.s the magnetic held is increased. The issue is 
that accurate numerical calculations become progressively more difficult the deeper one 
penetrates into the high-held/low-temperature regime, T^/|H| <C 1. (This is analogous 
to the difficulty of approaching extremality in the charged case, where the horizon 
temperature also vanishes.) By holding hxed El instead of Eb-, we are able to perform 
scans in B which avoid dipping too deeply into the very low temperature region. 

Another diherence concerns the dehnition of pressure (or stress) anisotropy. With 
our choice (2.27) for hxing the scheme dependent ambiguity in the stress-energy tensor, 
the resulting holographic relation (2.42), when evaluated at fi = 1/L, puts the trace 
anomaly solely in the transverse components of the stress, and So the pressure 
anisotropy (4.1), dehned as the diherence between transverse and longitudinal stress, 
when evaluated at /i = 1/L has a “kinematic” contribution of — plus a “dynamic” 
contribution of 3Kb4{v). In presenting results below, we will omit the uninteresting 
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Figure 14. Left: The subtracted/rescaled anisotropy function, Bs{v,u), as a fnnction of 
time V (in units of and inverse radial depth u, for the case of BL^ = 1.0, with initial 

Gaussian parameters chosen to match the initial pulse which generated fig. 5 (^ = 5 x 10“^, 
ro = 4, cj = 2 )) and energy density El = At late times the anisotropy function 

approaches the non-trivial profile of the eqnilibrium magnetic brane solntion discussed in 
section 2.7. Right: Corresponding evolution of the dynamical contribution to the relative 
pressure anisotropy , AV^yn/n^L, with both AV and el evaluated at the scale 1/L. The 
late time limit of the pressure anisotropy is non-zero, but too small to be easily visible, 
lim„^oo ARdyn(l/T)/KeL = 0.22. 


static kinematic contribution, and just plot the dynamic contribution 

APdyn = I + In 13\ (4.3) 

(relative to the energy density), evaluated at the renormalization point /r = 1/L?^ 

A further difference comes from the fact that equilibrium magnetic brane solu¬ 
tions are intrinsically anisotropic. The anisotropy function B does not vanish at late 
times, but rather settles down to the prohle of the equilibrium magnetic brane solu¬ 
tion discussed in section 2.7. This is illustrated in hgure 14, which shows the (sub¬ 
tracted/rescaled) anisotropy function Bs as a function of time v and inverse radial 
depth u. One sees similar features as in £g. 5: the initial pulse propagates outward, 
reflects off the boundary, and largely disappears into the horizon. But in addition one 
also sees that the anisotropy function is approaching the non-trivial profile of a static 
magnetic brane solution. 

^^Examining the holographic relation (2.42), one sees that simply shifting the renormalization point 
to /r = IL would accomplish the same removal of this uninteresting kinematic contribution to the 
pressure anisotropy. 
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Figure 15. Time dependence (in units of of the dynamical pressure anisotropy, at 

the scale of 1/T, for values of background magnetic field ranging from 0 to 4/L^. The energy 
density (at the curvature scale) is held fixed, and the parameters of the initial 

Gaussian pulse in the anisotropy are the same as in fig. 14. 

Correspondingly, the pressure anisotropy in the dual held theory asymptotically 
approaches a non-zero constant. To examine equilibration, it is the difference between 
the pressure anisotropy and its asymptotic value which is of interest. As a measure of 
(near) equilibration, the condition (4.2) is naturally replaced by 

[AV{t) - AV{oo)] /{ksl) < 0.1, (4.4) 


for all times t > t^q. 

The time dependence of the resulting (dynamical contribution to the) pressure 
anisotropy is shown in £g. 15 for a series of magnetic helds, BL^ = 0, 1, 2, 3, and 4. 
The energy density is held hxed ai El = \L~^ and the parameters of the Gaussian 
pulse in the initial anisotropy function are those of the pulse which generated £g. 5 
= 5 X 10 rg = 4, (T = |). For the hve cases shown, the ratios of magnetic held 
to the equilibrium temperature (squared) are given by B/T"^ = 0, 13.0, 30.2, 30.5, and 
26.3 for BL^ = 0, 1, 2, 3 and 4, respectively. Note that, at this hxed value of 
B/T"^ is not monotonic as a function of BL^. The energy densities at the intrinsic scale 
set by the magnetic held for this series of solutions are given by Eb/B’^ = cxd , 0.75, 0.36, 
0.36, and 0.39, respectively. 

Diherences in the late time values of the pressure anisotropy are obvious in hg. 15.^® 
^®The tiny positive late time pressure anisotropy barely visible in the BL^ = 1 curve is a consequence 


42 














bfP dyn/ 



Figure 16. Time dependence (in units of (vrT)”^) of the dynamical pressure anisotropy for 
values of background magnetic field ranging from 0 to 2.71/L^. The temperature is held fixed, 
ttTL = 1, and the parameters of the initial Gaussian pulse in the anisotropy are the same as 
in fig. 14. The values of for these solutions, in order of increasing field, are oo, 1.2, 

0.51, and 0.39, respectively. One sees very little sensitivity to the magnetic field in the time 
course of the response. 

Small temporal variations are also evident after v ~ 1.3. These are prodnced by the 
relaxation of the initial non-Ganssian backgronnd prohle of the anisotropy fnnction 
to the correct eqnilibrinm form. These small variations at relatively late times would 
be absent if we had constructed our initial anisotropy function by adding a Gaussian 
perturbation to the equilibrium solution (instead of merely adding a Gaussian to the 
leading asymptotic term). Given our choice of initial data, there are two distinct 
time scales in the equilibration shown in hgs. 14 and 15, the hrst associated with the 
boundary reflection and subsequent infall of the Gaussian pulse, and the second with 
the time it takes for the background anisotropy prohle to reach its equilibrium form. 
It is the latter which is responsible for the late time variations; fortunately, there is 
relatively little ambiguity in separating the two contributions to the dynamics. 

Once again, a notable feature of in the comparison of hg. 15 is the similarity in the 
time dependence of the pressure anisotropy during the pulse-driven period of signih- 
cant departure from equilibrium (0.2 < < 1.0), when the initial perturbation is 

rehecting from the boundary. Any reasonably dehned measure of equilibration time Gq 

of our removal of the static kinematic part of the anisotropy; the final value of the total pressure 
anisotropy, at the scale 1/L, monotonically decreases with increasing BL^ in this series of solutions. 
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of background magnetic field ranging from 0 to ^ for the same deep pulse (^ = 

ro = (T = used in fig. 9 and fig. 12. The energy density is held fixed. 

clearly does not vary much with magnetic held, and neither does tpeak for these rela¬ 
tively near boundary pulses. This insensitivity result relies, of course, on the constancy 
of the initial Gaussian perturbation in the anisotropy function, and also on our choice 
to hold hxed the energy density at the scale 1 /L. 

Fig. 16 compares the response, for different values of magnetic held, when one holds 
hxed the equilibrium temperature T (as well as the initial Gaussian perturbation), 
instead of hxing the energy density With time now plotted in units of (vrT)”^, one 
also sees remarkable similarity in the time dependence of the response, with the times 
of the hrst, second, or third peaks in the pressure anisotropy varying by less than 0.3% 
as B/T"^ varies from 0 to 26.7. 

Sensitivity to the magnetic held is signihcantly larger when the initial pulse is placed 
very close to the horizon. This is shown in hg. 17, which plots the time dependence of 
the dynamical pressure anisotropy for magnetic held values BL? = 0, 1, 1.5, 2, and 2.5, 
using the same “deep pulse” Gaussian parameters {A = = ^,<7 = ^), and hxed 

energy density = | which generated hg. 9. For this series of solutions we have 
\B\/T‘^ = 0, 13.0, 22.85, 30.16, and 31.95, and eB/B‘^ = oo, 0.75, 0.43, 0.36, and 0.35, 
respectively. 

With this deep initial pulse, diherences in the time course of the response are much 
more pronounced. Larger magnetic helds greatly suppress the size of the pulse-driven 
peaks in the pressure anisotropy (for a hxed amplitude of the initial Gaussian), and lead 
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Magnetic (S 7 ^ 0, p = 0) 


B/T'^ 

sb/T^ 

V/kT^ 


ReA/e^^ 



0 

73.06 

24.35 

0 

3.3521 ± 0.0001 

-2.9514 ±0.0001 

3.1195 

- 2.7466 i 

0.990 

72.98 

24.16 

-1.13 

3.357 

± 0.001 

-2.93 

±0.06 

3.124 

- 2.73* 

5.344 

80.74 

22.15 

-10.60 

3.372 

± 0.002 

-2.92 

±0.06 

3.217 

- 2.79* 

12.953 

125.85 

13.98 

-15.76 

3.264 

± 0.007 

-2.78 

± 0.01 

3.480 

- 2.96* 

17.821 

170.16 

3.79 

-9.36 

3.161 

± 0.002 

-2.69 

±0.04 

3.634 

- 3.09* 

22.836 

226.69 

-11.35 

5.00 

3.061 

± 0.008 

-2.60 

±0.03 

3.780 

-3.2H 

30.161 

328.21 

-42.21 

39.97 

2.94 

± 0.01 

-2.49 

±0.03 

3.98 

- 3.38* 


Table 2. Equilibrium energy density, average pressure, and pressure anisotropy, plus lowest 
quasinormal mode frequency, for various values of the external magnetic field. Reported 
values of energy densities and pressures are evaluated at a renormalization point [i = 
not at the (physically irrelevant) curvature scale. The pressure anisotropy is the complete 
valne, inclnding the — kinematic contribution which was removed in figs. 15 - 17. Results 
for the lowest quasinormal mode frequency are reported both in units of (middle colnmn, 
with uncertainties), and as well as in units of ttT (final column). 

to increasingly large and negative final values for the anisotropy. For the lowest curve 
with BL"^ = 2.5, the contribution of the Gaussian pulse is completely swamped by the 
contribution from the relaxation of the background prohle of the anisotropy function to 
the correct equilibrium form. Excluding this curve, the time tpeak of fhe hrst peak, as 
well the rough equilibration time tgq; characterizing the portion of the evolution arising 
from the Gaussian pulse, vary at most by 20% as the magnetic held ranges from 0 to 
2/L^. This is the largest difference in the relaxation time course we have seen in our 
exploration of magnetized plasmas with Gaussian initial perturbations. 

A constant temperature comparison (not plotted), analogous to hg. 16 but using 
the same deep pulse as in hgs. 9 and 17, shows variations in the time course of up to 
15% as \B\/T‘^ ranges from 0 to 22 — beyond which the response of the pulse cannot 
be clearly distinguished from the background evolution. 

To conclude this section, we report in table 2 equilibrium properties, plus our 
estimates for the lowest quasinormal mode frequency, for values of magnetic held which 
extend from small helds well into the strong held regime. The ratio B/T"^ ranges from 
zero to just over 30. The equilibrium energy density, pressure, and pressure anisotropy 
are given in units of T^, and have been converted from the fi = 1/L renormalization 
point used in our calculations (and above presentation) to the intrinsic scale jj, = \B\^/‘^. 
Results for the lowest quasinormal mode frequency are given both in units of and 
in units of ttT. These estimates are the result of htting the late time (4 < n < 9) 
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behavior of the pressure anisotropy to a decaying, oscillating exponential plus a constant 
equilibrium offset. To our knowledge, independent results from a linearized analysis 
of small fluctuations about the (numerically determined) magnetic brane solutions are 
not currently available. The leading quasinormal mode frequency varies by about 20% 
as B/T'^ ranges from zero up to 30, monotonically increasing with increasing held when 
measured in units of ttT, but slightly increasing and then decreasing when measured 
in units of . 

5 Discussion 

The above results show that to a good (often extremely good) level of accuracy: 

1 . the pressure anisotropy response is a linear functional of the initial anisotropy 
pulse prohle; 

2 . the time course of the response, measured in units set by the energy density,^^ 
is insensitive to the charge density or background magnetic held when the pulse 
prohle and the energy density are held hxed; 

3. the time course of the response, measured in units set by the equilibrium tem¬ 
perature, is insensitive to the charge density or background magnetic held when 
the pulse prohle and equilibrium temperature are held hxed. 

How can one synthesize these observations? Consider some feature in the time course 
of the response, such as the time of the hrst (or second) peak in the pressure anisotropy, 
or the approximate equilibration time discussed above. For simplicity, we will focus on 
the time tpeak of the hrst pressure anisotropy peak. This time must be some function 
of the equilibrium state parameters (energy density plus charge density or magnetic 
held), as well as the Gaussian pulse parameters (depth, width, and amplitude). 

Consider hrst the charged case. The response time fpeak is a function of the energy 
density e, the fraction x of the extremal charge density, and the Gaussian pulse param¬ 
eters tq, a, and A. But since the temperature decreases monotonically with increasing 
charge density (for hxed energy density) one may equally well regard the equilibrium 
state as labeled by the energy density e and temperature T, and write 


tpeak/L = f{eL'^,TL,ro/L,a/L,A) , (5.1) 

for some function / of the indicated arguments. We have written all quantities in 
dimensionless form using, in ehect, our computational units. Our results on the degree 

^^Specifically, £z, in the magnetic case. 
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of nonlinearity imply that the function / is nearly independent of the last argument, 
the Gaussian amplitude A. Only for our narrowest pulse, located right at the horizon, 
did the relative nonlinearity reach 10%. Away from this corner of parameter space, the 
nonlinearity was substantially smaller, rapidly falling to much less than a percent as 
the initial pulse becomes less localized at the horizon. So, to a good approximation, 
one can regard the function / as being independent of A. 

For narrow pulses, a tq, the dependence of the response time fpeak on the pulse 
width is negligible; there is a smooth cr —?■ 0 limit. To the extent that linearity is a good 
approximation, one may regard the response from wider pulses as superpositions of the 
response from narrower pulses (suitably arranged so that their sum reconstructs the 
wider pulse). As noted in the discussion of fig. 7, the first peak in the response is clearly 
associated with the propagation of the leading edge of the anisotropy pulse — the region 
of near-maximal slope on the side of the pulse closest to the boundary. Therefore, for 
pulses of non-negligible width a centered at some depth tq, one should expect that the 
time of the hrst peak in the response will be most similar to the corresponding response 
time for a narrower pulse located not at the depth tq, but rather at a depth of tq -|- na 
for some positive 0(1) multiplier n. At the same level of accuracy determined by the 
degree of nonlinearity, one should be able to merge the dependence of the response time 
^peak on the depth tq and width a of the initial pulse into a single effective depth rgg 
given by tq -|- na. (The accuracy of this simplification will be discussed below.) Hence, 
the above functional dependence for fpeak can be replaced by a simpler form, 

tpeak/L ^ g(eL^,TL,reff/L), (5.2) 

for some function g. Now, the results of section 4.2 (including hgs. 11 and 12 and 
associated discussion) show that the time course of the pressure anisotropy response has 
remarkably little dependence on the charge density when comparisons are made holding 
fixed the initial pulse parameters and the energy density. Since varying the charge 
density at fixed energy density is, as noted above, equivalent to varying the equilibrium 
temperature, this implies that the function g describing the response time (5.2) is 
nearly independent of the second argument, TL. At the same time, the comparisons 
at fixed temperature also discussed in section 4.2 (fig. 13 and associated discussion) 
show that the time course of the pressure anisotropy response also has remarkably little 
dependence on the charge density when the initial pulse parameters and the equilibrium 
temperature are held hxed. This implies that the function g describing the response 
time (5.2) is nearly independent of the hrst argument, Hence, at a level of accuracy 
determined by the minimal level of nonlinearity, and the minimal dependence on charge 
density in comparisons at both constant energy density and constant temperature, the 
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response time tpeak must be a function of only the effective depth of the initial pulse, 

^peak/-h ~ f , (5.3) 

for some function h. Finally, this function must be consistent with the scaling relations 
discussed in section 2.5, which imply that L^/r scales in the same fashion as a distance 
(or time) in the boundary theory. Consequently, the dependence of the response time 
on the effective depth must have the form 


tpeak ~ CL^/r^s , (5.4) 

for some dimensionless constant C. 

A similar line of reasoning is applicable to the magnetic case. Since we used = 
s{L) and BL^ as parameters labeling the equilibrium magnetic brane geometry in our 
comparisons of magnetic plasma response, it is convenient to view the response time 
^peak as depending on these parameters plus the Gaussian pulse parameters, 

fpeak/^ = BL'^,ro/L, ajL, A) , (5.5) 

for some function /. Once again, the observed near-linearity of the response allows us 
to simplify this to the form 


tpeak/L ^ g{eLL‘^,BL‘^,res:/L), (5.6) 

for some function g. The near-independence of the time course of the response on the 
magnetic held B, for hxed El and a hxed initial pulse, implies that the function g is 
nearly independent of its second argument. Because El does not transform homoge¬ 
neously under the scaling relations (2.43)-(2.45) (due to the use of the curvature scale 
L instead of a physical scale in the dual QFT for setting the renormalization point) 
consistency with the scaling relations requires that the function g be independent of 
its hrst argument and depend inversely on the third. So, just as for the charged case, 

tpeak ~ , (5.7) 

for some dimensionless constant C. 

Figure 18 compares this simple model with a sample of our data for neutral, 
charged, and magnetized plasmas. The left panel shows data for relatively narrow 
pulses of width a = 1/20, while the right panel shows data from rather wide pulses with 
width cr = 1/2. The abscissa for both panels is the inverse effective depth Ues = l/^'efr- 
(Data points at the largest values of Ues shown in these plots come from pulses centered 
very near the horizon.) In both panels, rather good agreement with the simple model 
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Figure 18. Time of the first peak in the pressure anisotropy response as a function of the 
inverse effective depth of the pulse, Ues = l/reff- The left panel shows results for narrow 
pulses with a = 1/20, while the right panel shows results from wide pulses with a = 1/2. 
Blue circles represent data points from neutral (uncharged, unmagnetized) plasmas, maroon 
squares represent data points from charged plasmas at 80% of extremality, and gold diamonds 
represent data points from magnetized plasmas with BL^ = 1.5. The green straight line shows 
the prediction of our simple model (5.7) with n = 2.5 and C = 2.04. 


is found when the multiplier n in the effective depth Ves = + na is chosen to be 2.5, 

and the coefficient 

C^2M. (5.8) 

For both sets of data, the model is least accurate for pulses located very close to the 
horizon. Accuracy for the case of magnetized plasma is a bit worse than for charged or 
neutral plasma. But for all cases, even in the near-horizon regime, this simple model 
works at about the 20 % level or better. As pulses move away from the horizon, the 
accuracy rapidly improves. 


6 Conclusion 

In weakly coupled plasmas, adding a conserved charge density to the system (e.g., flavor 
charge in a QCD plasma) signihcantly changes screening lengths. When the associated 
chemical potential is of order vrT, the relative change in the Debye screening length is 
0(1) [51]. Such changes in screening lengths signihcantly affect transport coefficients 
like the viscosity [52]. Lattice studies [53] of the effect of a baryon chemical potential in 
the deconhned phase, when the system is not asymptotically weakly coupled, claim to 
hnd measurable sensitivity in the Debye screening length, comparable to perturbative 
estimates, but with quite large error bars. Lattice QCD studies of magnetoresponse 
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[29-32] also find substantial changes in thermodynamics when the magnetic held energy 
density becomes large compared to T^. 

Consequently, when this work on strongly coupled A/” = 4 SYM plasma was initi¬ 
ated, we expected to hnd signihcant changes in equilibration dynamics when a conserved 
charge density is added to the plasma, or when the system is placed in a background 
magnetic held. The most notable result we have found is that this expectation was 
wrong. At least within the range of charge densities we studied, up to 80% of ex- 
tremality, the equilibration dynamics is remarkably insensitive to the presence of a 
conserved charge density. Additionally, magnetic helds which are well into the strong 
held region, B/T'^ ^ 1, induce almost no change in the equilibration time course. 

Ehorts to use results of holographic calculations in strongly coupled A/^ = 4 SYM as 
the basis for predictions about real heavy ion collisions [21] are inevitably hampered by 
our limited understanding of the ehect on the relevant dynamics of changing the theory 
from real QCD to a supersymmetric Yang-Mills model theory. From this perspective, 
the insensitivity of the equilibration dynamics to the charge density is reassuring, as 
this provides an example where changes in the plasma constituents have very little 
impact on the overall dynamics. 

A further notable feature in our results is the remarkably small degree of nonlin¬ 
earity in the dynamics governing the pressure anisotropy. Despite the fact that one is 
solving the highly nonlinear Einstein equations, the dependence of the induced pres¬ 
sure anisotropy on the initial anisotropy function is surprisingly close to linear. We 
hnd deviations from linearity of at most ~ 30%, and this only for initial disturbances 
which are crafted to reside very deep in the bulk. (This is consistent with earlier work 
in ref. [9].) For disturbances localized even modestly above the horizon, the degree of 
nonlinearity quickly drops to sub-percent levels. This near-linearity holds even when 
the system is far from equilibrium, with pressure anisotropies which are large compared 
to the energy density. 

In the charged case, the sensitivity of the lowest quasinormal mode frequency to the 
charge density is signihcantly larger than the sensitivity we hnd in, and shortly after, 
the far-from-equilibrium period of the equilibration dynamics. The lowest quasinormal 
mode dominates the equilibration process at sufficiently late times when all higher 
modes have decayed away and become negligible relative to the lowest mode. 

The near-linearity which we hnd in the dynamics implies that the deviation from 
equilibrium, even during the far-from-equilibrium portion of the process, can be rep¬ 
resented quite accurately as a sum of quasinormal modes obeying linearized small 
huctuation equations. The minimal sensitivity of the equilibration dynamics to the 
charge density, substantially less than the sensitivity of the lowest quasinormal mode 
frequency, suggests that during most of the equilibration process many quasinormal 
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modes above the lowest one are contributing, with the sensitivity to the charge density 
quickly falling with increasing mode number. This is something which could be tested 
directly in a linearized analysis but, to our knowledge, has not yet been done. 

In the magnetic case, we hnd only modest (few percent) changes induced by the 
magnetic held in the time course to equilibration, and in the lowest quasinormal mode 
frequency characterizing very late time behavior. It would be nice to have independent 
calculations of quasinormal mode frequencies for magnetic branes, including studies of 
the held dependence of higher mode frequencies. 

Although we no longer have reason to expect large ehects, it would be natural to 
generalize the study of homogeneous equilibration to the case of plasmas with both 
non-zero charge density and a background magnetic held. Extensions to more com¬ 
plicated inhomogeneous settings, such as colliding shock waves, are also of interest. 
It would be desirable to gain a clearer understanding of the connection between our 
choices of gravitational initial data on null slices and boundary observables such as 
multi-point stress-energy correlators. Can more operational procedures, such as time- 
dependent background helds [7], produce far-from-equilibrium states which resemble 
those produced by our “deep pulse” initial data? We hope future work can shed light 
on some of these topics. 
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